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C/3 , ABSTRACT 

,'^, - A rapid timing analysis of VLT/ULTRACAM (optical) and RXTE (X-ray) observations of 

the Galactic black hole binary GX 339-4 in the low/hard, post-outburst state of June 2007 

is presented. The optical light curves in the r', g' and u' filters show slow (^ 20 s) quasi- 

^ ' periodic variability. Upon this is superposed fast flaring activity on times approaching the 

lO ! best time resolution probed (~ 50 ms in r' and g') and with maximum strengths of more 

OO . than twice the local mean. Power spectral analysis over ~ 0.004-10 Hz is presented, and 

>0 ' shows that although the average optical variability amplitude is lower than that in X-rays, the 

peak variability power emerges at a higher Fourier frequency in the optical. Energetically, we 
If^ • measure a large optical vs. X-ray flux ratio, higher than that seen on previous occasions when 

^^ ' the source was fully jet-dominated. Such a large ratio cannot be easily explained with a disc 

^D . alone. Studying the optical:X-ray cross-spectrum in Fourier space shows a markedly different 

behaviour above and below ^^0.2 Hz. The peak of the coherence function above this threshold 
is associated with a short optical time lag with respect to X-rays, also seen as the dominant 
feature in the time-domain cross-correlation at «150 ms. The rms energy spectrum of these 
/\ ' fast variations is best described by distinct physical components over the optical and X-ray 

H I regimes, and also suggests a maximal irradiated disc fraction of 20 per cent around 5000 A. 

If the constant time delay is due to propagation of fluctuations to (or within) the jet, this is 
the clearest optical evidence to date of the location of this component. The low-frequency 
quasi-periodic oscillation is seen in the optical but not in X-rays, and is associated with a low 
coherence. Evidence of reprocessing emerges at the lowest Fourier frequencies, with optical 
lags at ^10 s and strong coherence in the blue u' filter. Consistent with this, simultaneous 
optical spectroscopy also shows the Bowen fluorescence blend, though its emission location 
is unclear But canonical disc reprocessing cannot dominate the optical power easily, nor 
explain the fast variability. 

Key words: accretion: stars - individual: GX339-4 - stars: X-rays: binaries. 



1 INTRODUCTION of time scales (e.g. Ivan der Klislll989l ). This is a manifestation of 

rich underlying structure in the accretion environments of these 
The radiative output of active X-ray binaries (XRBs) peaks at X- sources on a range of physical scales. X-ray timing studies can 

ray wavelengths, hence their name. A characteristic property of thus provide key information for our understanding of volatile 

this radiation is its strong and aperiodic variability over a range 
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XRB environme nts, and X-ray timing is now a mature field in 
this respect (e.g. Hasinger & van der Klisll989l : lvan der Klisl 19951 : 



iRemillard & McClintockll2006h 

At lower energies, e.g. in the optical and infrared, the picture 
is thought to be complex. For many low mass XRBs, the donor 
star contribution can be ignored during active accretion states. Op- 
tical emission from the accretion disc arises predominantly from its 
outer, cool portions with temperatures much below ~10® K. Con- 
sequently, large variations of this component should occur only on 
a relatively- slow viscous timescale. Any faster variations may be 
the result of high-energy photons being reprocessed on the outer 
disc. In fact, such X-ray heating is thought to be the primary en- 
gine generat ing the bulk of the ultraviolet and o ptical emission seen 
from XRBs jvan Paradiis & McClintoc3ll994) . Weaker, but highly 
significant reprocessing on the companion star is also an impor- 
tant effect which can be used to determine the physical parame- 
ters of the biriary sy stem (e.g. iHomdl 19851 : iMarsh & Homelll988l : 
lO'Brien et ai]|2002l) . 

Observations are now also beginning to provide evidence 
of fast variability (on ~1 s timescales or less) at low energies 
which is not associated with reprocessing, and recent works at- 
tribute much of this to non-thermal contributions from the jet and 
corona. Some observations of fast low energy varia tions in sev- 
eral sources have actually existed for many ye ars (e.g. lMotch et al.l 
ll982l : IUemura et al.ll2002l : lFender et al.lll997h . In several cases, si- 
multaneous low and high energy observations argue against a re- 
processing origin. This was best demonstrated in XTE Jl 1 18-1-480, 
for which rapid optical and X-ray timing observations revealed 
complex flux variations, including an anti-correlation between op- 
tical and X-rays - exactly the opposite of that predicted by repro - 
cessing scenarios i Kanbach et alj|200l1 : ISpruit & Kanbachll2002h . 
Many works have successfully modelled the optical power and 
broad-band spectrum as a result of (cyclo-)synchrotron emission, 
though the location and char acteristics of the emitting plasma re- 
main a matter of debate (e.g. Merloni et alj200CI : lEsin et alj200ll : 
iMarkoff et alj2005l : lYuan et alj200^ . and more). Self-consistently 
incorporating the timing properties in this framework has proven 
to be trickier. An energy reservo ir model feed i ng a common jet 
and corona system developed by Malzac et al.l j2004h can repro- 
duce both the shape and the strength of the observed timing corre- 
lations very well. But physical modeling of the reservoir itself and 
of the impulses of energy release governing the observed variability 
remains somewhat ad hoc, being modell ed simply as a linear super- 
position of fluctuating stochastic shots. iGandhil j2009t) discusses a 
simple modification of this phenomenological model to reproduce 
another observable (the 'rms-flux' relation), by making the light 
curves a non-linear coupling of shots, rather than an additive super- 
position. 

Other sources with intriguing multi- wavelength timing corre- 
latio ns include GRS 19 15-1-105, studied bv lEikenberrv et ail J1998l) 
and I Aral et alj 120091) . among others. And in a recent study , 
lOandhi et alj |2008| . hereafter Paper I) and lOurant et alj ^200^ 
observed two other Galactic BH candidates in the X-ray low/hard 
state, GX 339-4 and SWIFT J1753.5-0127, simultaneously in X- 
rays with the Rossi X-ray Timing Explorer (RXTE) and in the 
optical with ULTRACAM mounted on the Very Large Telescope 
(VLT). Cross-correlating the light curves in the two waveband s re- 
vealed complex variability patterns in both systems. Dura nt et al.l 
( 120091) have discussed the second source in detail, suggesting that 
optical cyclotron anti-correlates with an X-ray (presumably in- 
verse Compton dominated) response (e.g. iFabian et alJI 19821) . In 



the present papeo we present the full spectral and timing prop- 
erties of GX 339-4, supplementing the cross-correlation function 
presented in Paper I. 

GX 339-4 (4U 1658-48) is one of the best-studied Galac- 
tic black hole (BH) candidates over a wide range of wavelengths. 
The source is a classic X-ray transient. It is usually found in 
the active, low/hard flux state, but is known to undergo dra- 
matic changes on month-year timespans, showing a broad range 
of X-ray states from 'qui escence', to the 'very high' state (e.g. 
Izdziarski et alj|2004 : Dunn et alj|2008i) . S uch changes are ac com- 
panied by strong modulations i n its radio JCorbel et alJuOOOf) . and 
optical JMakishima et al.lll986f) properties. The source has proven 
to be key for the discovery and elucidation of several important 
aspects of XRB behaviour. It was one of the first to show clear 
evidence of the existence of an in ner accretion disc in the X-ray 
low/hard state JMiller et alJuOOq) . It is the prototype source for 
defining an important relationship between the radio and X-ray 
fluxes of XRBs (e.g. Corbel et al. 2000: Gallo et al. 2003). 

It was also o ne of the first sou rces found to show very fast 
optical flickering JMotch et al.lll982h . Follow-up rapid simultane- 
ous optical and X-ray studies uncovered evidence of a slow anti- 
correlated response between the two bands, but such studies re- 
mained scarce (Motchetal. 1983: Makishima et al. 1986). In Pa- 
per I, we presented a new correlation feature in the faint low/hard 
state, with a sharp, sub-second optical response in a red (r') filter, 
delayed with respect to X-rays. As was discussed in that paper, and 
elaborated upon herein, simple linear reprocessing transfer func- 
tions cannot explain this feature. We now present our full multi- 
filter variability dataset, and relate this with the simultaneous X-ray 
timing properties over a range of timescales. The optical and X-ray 
power spectra and correlations in time and Fourier space are anal- 
ysed. Supplementary optical spectroscopy with the FORS2/VLT in- 
strument is also presented. The observations are discussed in the 
context of multiple varying components with some (as yet not fully 
understood) underlying interactions. Our broad-band analysis pro- 
vides first simultaneous constraints on the physical characteristics, 
emitted powers and relative locations of these components within 
the accretion environment. 

The donor star of this low mass XRB system is very faint 
JShahbaz et alj2001r) . which complicates measurement of physical 
parameters of the binary, but means that the optical flux is com- 
pletely dominated by the accretion flow, which is our main inter- 
est herein. We assume a compact object mass of 6 Mq, consis- 
tent with the best determ ination of the mass function of 5.8±0.5 
M^ jHynes et al.i 2003b), and a distance of 8 kpc (a range of 6- 
15 kpc has been sugg ested in the literature: e.g. Hynes et al. 2004|; 
Izdziarski et al.ll2004h . The most-widely accepted value of the bi- 
nary orbital period is 1.75 days, based on outburst radial velocity 
variations of th e optical Bowen em ission blend reprocessed on the 
secondary star l IHvnes et alj2003bh . and we adopt this value here. 



2 OBSERVATIONS AND DATA ANALYSIS 

2.1 Optical timing (ULTRACAM) 

Our optical observations were carried out with ULTRACAM, 
which is a fast ima ging camera employing frame-transfer CCDs 
I IDhillon et al.l2007h . Although not a 'photon-counting' instrument. 



^ Ba.sed on observations obtained during ESO programs 079.D-0535 and 
279.D-5021. 
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ULTRACAM satisfies two critical requirements for astronomical 
studies with CCD detectors on fast (sub-second) timescales - 1) a 
low readout noise and, 2) a small dead time between exposures. Un- 
like proportional counters, CCD dead time does not depend on the 
incident count rate. Instead, it is just the time taken to read out, or 
clock, the exposed CCD pixels, during which any incident photons 
are not accumulated. Frame-transfer CCDs are constructed with an 
extra storage area identical to their exposed pixel area, and clock- 
ing quickly transfers the photo-generated charge to this storage area 
from where it is subsequently read out and digitised. The maximum 
dead time in ULTRACAM is only 24 ms. But this holds true for full 
CCD read out and signal digitisation. The dead time can be much 
reduced by reading out only small portions of the CCD which are 
of interest. 

UUTRACAM was mounted on the Nasmyth focus of 
VLTAJT3 as a visitor instrument during 2007 Jun 9~24, and our 
observational time window was governed by this mounting period. 
We happened to catch GX 339^ in its typical X-ray low/hard 
state, soo n after it had emerg ed from an outburst that peaked in 
mid Feb ( iKalemci et alj 120071) . Fig. [T] shows the long-term light 
curve of GX 339-4; our days of observations are denoted by ar- 
rows. We chose to monitor the source for four hours of integra- 
tion, split equally across four nights. No observing was possible 
on the night of UT Jun 12 (Chilean night starting Jun 11) because 
of clouds. Weather steadily improved thereafter, with changeable 
conditions on Jun 14 and 16, and very stable conditions on Jun 18. 
These nights are hereafter referred to as Nights 4, 3, 2 and 1, re- 
spectively, ranked according to the weather. 

ULTRACAM has a 2.6x2.6 arcminute field at the VLT. In or- 
der to achieve a reasonable signalmoise, the timing resolution of 
our observations was decided interactively at the beginning of each 
sequence of exposures based on real-time weather conditions. The 
size of the CCD window to be read out was adjusted in order to 
achieve this time resolution, resulting in small window sizes on 
all nights (not larger than 100x100 pixels centred on the target). 
In addition, to achieve the fastest frame rates, we used the 'drift' 
observing mode. Instead of clocking (i.e. charge transferring) the 
entire exposed area of the CCD to the unexposed area, drift mode 
clocks only the pixel rows corresponding to the window region. 
As clocking is typically much faster than signal digitisation, the 
dead time (equal to th e clocking tim e) is very small. The interested 
reader should refer to lDhillon et al.l 12007) for details. For our ob- 
servations, the frame rates ranged from just under 50 ms to about 
136 ms, and dead times were approximately 1.6, 2.3 and 2.3 ms on 
Nights I, 2 and 3, respectively. These correspond to about 3.2, 1.7 
and 1.7 per cent of the cycle time on each night. Table [T] lists the 
observation log. Note that the cycle and frame times stated therein 
are average values and the standard deviation of these is ~ 5 /is, 
which means that the light curves are not strictly regular at the sub- 
millisecond level. But this is not a problem because of the accurate 
time sampling. 

Source flux was monitored by observing a bright comparison 
star (USNO-Bl.O 0412-0564199) within the field of view, centred 
on a second CCD window identical in size to that used for the tar- 
get. ULTRACAM is equipped with three separate, but identical, 
detectors, and with the use of beam-splitters, fully simultaneous 
observations in three separate optical filters can be carried out. We 
used the r', g' and u SDSS filter set. Note that low atmospheric 
(and inter-stellar) transparency in u necessitates co-addition of 
many frames for source detection, with the result that cycle times 
of faster than a few seconds were not probed for this filter. 




200 



Figure 1. Long-term RXTE All Sky Monitor light curve of GX 339- 
4, obtained and rebinned from the quick-look results provided by the 
ASM/RXTEteam (http://xte.mit.edu/asmlc/). The arrows are the nights of 
our observations, denoted Nights 4, 3, 2 and 1 from left to right, in terms 
of improving weather. The first day of several months in 2007 are marked 
at the top. As is appaient, our observations were caiTied out soon after the 
source had returned to a low flux (but active) state following its 2007 out- 
burst. 



The ULTRACAM pipelinqj was used for data cube reduction, 
including bias subtraction, flat-fielding and photometry, for each 
of the three filters separately. Aperture photometry was carried out 
by profile fitting and optimal extraction of target and comparison 
counts separately, and annular background subtraction. All source 
light curves were divided by the corresponding comparison star 
light curves normalised to their mean. The photometric solution 
is very stable on Night I, but degrades progressively on the other 
two nights, with the source being lost on several frames during our 
hour-long monitoring on Night 3. The pipeline provides some con- 
trol by allowing variable aperture sizes and centering options, but 
the worst of these bad weather periods - a single contiguous 25 min 
segment on Night 3 (see Fig.|2j - has been avoided. Relative pho- 
tometry was thus important for Nights 2 and 3, especially for the 
long-term (> minute timespan) variations. Additionally, there are 
intermittent source detection drops for the lower sensitivity u data 
throughout both these nights, and hence these data are not analysed. 
For Night 1, 50 u' co-adds were used, resulting in best time bin of 
~2.5 s for this filter. 

In short, we use the full r' and g' datasets on all three nights 
for the timing analysis presented in this paper, except for the bad 
transparency segment on Night 3. For the u' filter, only data from 
Night 1 are used. 



2.2 Optical spectroscopy (FORS) 

We were also awarded director's discretionary time for short 
optical spectroscopic observ ations with the F0RS2 instru- 
ment ( App enzeller et alj|l998h . simultaneous with ULTRACAM. 
Changeable weather conditions meant that only two spectra were 
obtained simultaneously, with the data from the best night (Night 
1) presented herein. The GRIS 300V grism with the GG375 or- 



http://deneb.astro.warwick.ac.uk/phsaap/software/ultracam/html/ 
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Figure 2. Times of the optical and X-ray observations. For botli telescopes, 
the marked intervals denote the effective 'good' time intervals only. The 
ULTRACAM observations have small (few second gaps) which are ex- 
cluded on Nights 1 and 2 from the analysis. The worst period of highly 
variable transpai'ency during Night 3 is a hatched 25 m segment starting at 
UT06:08, which was also excluded. 



der sorting filteu was used for a total integration time of 200 s, 
starting at UT07:09. The effective wavelength coverage was over 
the ~ 4000-7500 A regime, with a resolving power i?~700 at 
5000 A. Standard data reduction, wavelength calibration and spec- 
trum extraction were carried out in IRAF. Flux calibration was car- 
ried out against a spectro-photometric standard star (LTT 7987; 
iHamuv et alJI 19941) observation from the same night. 



2.3 X-rays (RXTE) 

In X-rays, we observed GX 33 9-4 for four hours (Observation ID 
prefix 93119-01) with RXTE JBradt et al.|[l993h . coordinated to 
be closely simultaneous to the optical observations over the four 
nights. In practice, RXTE orbital constraints resulted in the exact 
length of overlap being slightly shorter than the desired 60 mins on 
each night. Table[T]lists the RXTE observing log and Fig.[2]shows 
the final overlap period of the optical and X-ray observations in 
graphical form. 

RXTE has two narrow-field instruments. The one with a low 
energy (~2-6 keV) response is called PCA (the Proportional 
Counter Arrav; |jahoda et alJll996r) . and that sensitive at higher en- 
ergies (~ 15- 250 keV) is HEXTE (t he High Energy X-ray Timing 
Experiment; iRothschild et al.ll 19980 . HEXTE employs two detec- 
tion clusters, each originally capable of rocking between the source 
and a background position. At the time of our observations, one of 
the two (cluster A) could no longer operate in the rocking mode and 
hence real-time background subtraction was not possible for this. 
Furthermore, the count rate of GX 339-4 in HEXTE was very low 
and dominated by the background; the net rate was only 3 ct s^^. 
Hence we do not analyse the HEXTE timing data further herein. 
But the net counts from cluster B accumulated over the full ex- 
posure duration were sufficient for extraction of a spectrum with 
useful signal. 



■^ FORS2 User's Manual (http://www.eso.org/sci/faciHties/paranal/instru 
ments/fors/doc/) 



The PCA has five proportional counter units (PCUs), and two 
or three of these were operational during our observations, the 
number changing from night to night. For data reduction, we fol- 
lowed standard procedures described in the RXTE cookbook and 
other material provided by the J?XTE Guest Observer Facility, us- 
ing FTOOL^ telackbumlll995h and HEASOFT v.6.8. The detected 
photons are analysed and delivered in several formats. In the event 
mode recommended for weak sources, GoodXenon mode, every 
event is telemetered to the ground (with intrinsic time resolution of 
~I fis) and this was the mode used for our timing analysis. Spectra 
were extracted from both instruments. For the PCA, we used the 
binned Standard2 mode, with a time resolution of 16 s. Basic 
good time intervals were selected by filtering on a minimum ele- 
vation angle of 10 degrees from the Earth's limb and a maximum 
pointing offset of 0.02. Each PCU has three Xenon detector layers, 
all of which were summed. Background files were generated by 
running the pcabackest utility, which estimates the background 
using Standard2 data files and synthetic background models. 
For HEXTE, data and background were extracted from the good 
channels of cluster B. Standard FTOOLS routines were used for re- 
sponse matrix generation and dead-time correction. Spectral mod- 
elling was carried out over ~3-30 keV (PCA) and 20-200 keV 
(HEXTE). Systematic uncertainties of 0.5 per cenjfland I per cenjj 
were included for the two instruments, as recommended. 



2.4 Barycenterlng of light curves 

Several of the timing results that we will describe below depend on 
accurate calibration of any time lag between the two light curves. 
Hence, transformation of the light curves to a common time frame 
is an important part of the data reduction process. In X-rays, this 
was accomplished with the f axbary routine to convert the RXTE 
data file times from mission time to a Bary centric frame in Ter- 
restrial Time (TT). In the optical, a custom-built routine called 
tcorr from the ULTRACAM pipeline was used, which carries 
out absolute time conversions using the SLA C library for abso- 
lute timing (Pat Wallace, priv. comm.), with inputs being the source 
coordinates (we used 12000 epoch coordinates as RA=17:02:49.5 
and Dec=-48:47:23), the name of the telescope (whose latitude, 
longitude and altitude are internally pre-defined for the VLT as 
70°24'9".9W, 24°37'30".3S and 2635m, respectively) and the ob- 
serving time in UT for all time bins. Absolute timing in ULTRA- 
CAM is achieved by GPS time s tamps from multiple satellites, and 
is good to at least 1 ms (.Dhillon et al.ll2007l) . In the fast drift ob- 
serving mode, there are additional corrections required for time 
stamping, because an exposed CCD window is not read out (i.e. 
digitised) immediately upon completion of its exposure. Instead, 
with each exposure, it is clocked by the number of rows it contains 
and is read out only when it reaches the CCD edge. But it is a sim- 
ple matter for the pipeline to determine the time stamp based on the 
(very accurately known) clocking and digitisation times per pixel. 
In X-rays, the absolute accuracy of the raw times are better than 
100 /lO The accuracy of the Barycentre transformation codes has 
been compared against pulsar timing, and found to be accurate to 
better than 100 /is. Similar transformation applied to the data for 



* http://heasarc.gsfc.nasa.gov/ftools/ 

^ http://www.universe.nasa.gov/xrays/programs/rxte/pca/doc/rmf/pcarmf- 

11.7 

^ http://heasarc.gsfc.nasa.gov/docs/xte/ftools/xtefaq_answers.html 

^ http://heasarc.gsfc.nasa.gov/docs/xte/abc/time.html 
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anothe r source - SWIFT J1753. 5-0127, described in lDurant et alj 
( 120081) - gave an optical vs. X-ray cross-correlation with a sharp 
change exactly at zero delay, giving extra confidence in the timing 
calibrations. 



3 RESULTS: SPECTROSCOPY 



3.1 Optical 



The observed optical spectrum is shown in Fig. [S] GX 339-4 
lies behind large interstellar reddening, whose exact value has 
been debated in the literature (see, e.g., extensive discussion in 
Izdziarski et alj Il998r) . Among the lat est works, [Z dziarski et al. 
( Il998l) report gR-v=1.2±0.1. while iBuxton & Vennesl J2003I) 
quote a value of ~ 1.1 ±0.2. Assuming the standard Galactic value 
of the totahselective extinction _Rv=3.1 gives ^v~3.5 mags, 
which was adopted for dereddening corre ctions based on the Galac- 
tic extinction law of lCardelli et alj ( 119891) . The error on the redden- 
ing translates into a relatively-large uncertainty on the correspond- 
ing extinction Aylv of at least 0.5 mag, and hence a dereddened 
l/-band flux uncertainty of 1.6, which is likely to be a lower-limit 
because _Rv is also known to vary along different sightlines. The 
final dereddened and flux-calibrated spectrum is also shown in the 
figure. 

The observed spectrum possesses a continuum that rises to- 
wards red wavelengths faster than a flat spectrum source, over- 
laid with several emission lines. The dereddened continuum, on the 
other hand, rises towards blue wavelengths, with a simple power- 
law parametrisation having a form _Fa~A~^'*, in units of flux den- 
sity per unit wavelength. A systematic uncertainty of about 0.7 in 
slope is possible for the full range of reddening correction uncer- 
tainty above. But, as can be discerned from the figure, a single 
power-law is not a good fit, with the slope below ~5000 A being 
bluer (Fa~A~^°), and that above being redder (Fa~A~^°) than 
expected from a single power law. Such a 'convex-shaped', broken 
power-law continuum is not a result of over-correcting for the red- 
dening; we confirmed this by using Av values at the lower end of 
expected range. We also note that the flux calibration fit resulted 
in an excellent rms of 0.01 mags over the range of 4000-7500 A, 
and a broken slope appeared on the other night of data as well (not 
shown). 

The observed spectral energy density normalisation is 
Fa (5500 A)= 6x10^^** erg s"^ cm"^ A"\ which is equivalent to 
Vvcga ~ 17. This matches contemporaneo us photometric monitor- 
ing of the source ( IBuxton & Bailvnll2007r) . Our average extinction 
correction and related systematic uncertainties then yield Mv=- 
1.02±0.5 or ALa(5500 A) « 6+2 x 10^* erg s"\ for the assumed 
distance of 8 kpc. 

The strongest observed emission lines include the Balmer se- 
ries from Ha to H5, the Bowen blend around 4640 A and several 
He I and He II lines. The se lines are typical for GX 339^ in the 
low/hard state dSoria et al.|[l999l : IBuxton & Vennesll2003[) . See Ta- 
ble[2]for the line strength measurements. Most lines, especially the 
strong ones (e.g. Ha and He I A4686) appear single-peaked within 
the limits of our modest spectral resolution. The Na I interstellar 
absorption line is also visible for which we measure an equivalent 
width of 3.8 A, though no correction for optical depth effects has 
been made. This lies with in the range of line strengths tabulated by 
IBuxton & Vennesl ( l2003h . 



Line 


Eq. Width 

A 


HS 4102 
H7 4342 

Bowen blend (BB) 
He II 4686 


1.00±0.26 
1.09±0.24 
2.12±0.39 
4.61±0.52 


H/3 4861 
He I 4922 
He I 5016 
He II 5411 


1.86±0.21 
0.56±0.08 
1.23±0.24 
0.70±0.35 


He I 5876 
Na I 5891 


0.90±0.22 
-3.84it0.24 


Ha 6563 
He I 6678 


6.17±0.17 
0.75±0.17 



Table 2. Line features detected in the FORS optical spectrum of GX 339^ 
obtained on Night 1. Eq. Width is the observed equivalent width measured 
from a single symmetric gaussian fit, with negative values denoting absorp- 
tion. Error propagation on the line width, normalisation and continuum level 
was used to determine the stated Icr uncertainties. 



3.2 X-ray 

X-ray spectroscopy quasi-simultaneous with our observations has 
been discussed in detail by Tomsi ck et al. (2008) - their 'spectrum 
2' was last observed on 2007 Jun 14, corresponding to our Night 
3. They used data from the Swift satellite and its X-Ray Telescope 
(XRT) instrument, in combination with RXTE. The XRT employs 
a CCD detector and good energy resolution below 10 keV, making 
it superior to RXTE/PCA in terms of isolating sources and for mod- 
elling the low energy (and Fe Ka) spectrum regions. So we present 
only minimal details of our spectroscopic analysis here. 

A simple absorbed power law with a photon index 
r=1.63±0.01 (uncertainty for Ax^=l) absorbe d by a fixed Galac- 
tic iVH=5.3xlO^^ cm~^ dOickev & LockmaiJlT99a) gives a best 
fit statistic of x'^/dof=99.5/96 (dof being the number of de- 
grees of freedom) over 3-200 keV, and no significant cross- 
calibration offset between the two instruments. Adopting a re- 
flection model (PEXRAV) instead, results in moderate improve- 
ment at x^/dof=96.1/95, for r=1.66±0.04, reflection fraction 
f2/27r=0.15±0. 1 but with unco nstrained power-law c ut-off energy. 
This is in good agreement with Tomsick et al.l ( l2008h . and we refer 
the reader to that paper for further discussion on the spectral mod- 
elling. Here, we simply measure the source power from our data. 

The observed source 3-10 keV flux from the power law model 
is _Fli_io=1.39(±0.01)xlO"^° erg s'^ cm~^ where the error is 
for 90 per cent confidence, estimated using the cf lux model pa- 
rameter in XSPEC using the PCA data alone. Extrapolating the 
model energy range and correcting for absorption gives a 2-10 
keV flux of 1.75(±0.02)xl0"^° erg s"^ cm"^ and a luminosity 
of L2-io=1.3x 10^® erg s~^ at 8 kpc. For comparison with the 
Eddington luminosity (LEdd), a broader band is more appropri- 
ate; our model gives Li_ioo=5.3x 10"^*' erg s^^, corresponding to 
0.007 X I/Edd. No important luminosity variations are seen between 
the different nights. 



3.3 Broad-band spectral energy distribution 

Fig. m shows the broad-band spectral energy distribution (SED) of 
GX 339^ including the simultaneous X-ray and optical data from 
the previous section. Night 1 measurements are used here, as this 
had photometric weather conditions. We note that the optical spec- 
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Night 


Date 
2007 Jun 


RXTE 






ULTRACAM 




Overlap 




Times 


Active PCUs 


Times 


Exp time 


Cycle time 






UT 


UT 




IDs 


UT 


ms 


ms 


min 


4 


12 


04:51-05:54 




024 




_ 






3 


14 


05:31-06:37 




02 


05:43-06:59 


134.0 


136.3 


49 


2 


16 


04:40-05:50 




024 


04:54-06:22 


131.0 


133.3 


47 


1 


18 


06:55-08:03 




02 


07:03-08:00 


48.0 


49.6 


56 



Table 1. Observation log. The first column states our night designation, based on weather conditions. The third column lists the RXTE good time intervals. 
Night 4 was completely clouded out, and is not discussed in this paper. The 'Exp time' and 'Cycle time' columns list the individual frame exposure times, 
and time difference between successive frames, respectively, rounded off to tenths of a milli-second. For the u' filter, only the data from Night 1 are analysed, 
and these required 50 frame co-adds, resulting in correspondingly longer best time resolution. The last column states the length of the overlapping optical and 
X-ray observational interval, rounded off in minutes. There are a handful of short gaps in the optical and X-ray data on different nights (and a longer period of 
bad transparency on Night 3; see Fig.|2), which were removed from the analysis. Only contiguous sections of specific length (typically 256 s) were used for 
power spectra and coherence computations, as desciibed in the text, leading to further slight shortening of the overlap time in these cases. 
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Figure 3. Optical spectrum of GX 339^ obtained with VLT/F0RS2 on Night 1. (Left) The observed flux-calibrated spectrum, with main line identifications 
labelled. The circled cross denotes telluric absorption. (Right) The dereddened spectrum (assuming Av=3.5) plotted in XFx units and against frequency for 
better comparison with the broad-band spectral energy distribution discussed later. Note the lower limit of the y-axis is above zero. The dotted lines indicate 
the continuum power law fits below (blue) and above (orange) 5000 A, extrapolated over the full frequency range. 



tral fluxes are close to those obtained by conversion of the mean 
count rates in ULTRACAM monitoring to BV filter fluxes, assum- 
ing standard zero-points and fi lter transformatio ns applicable for 
the ULTRACAM Sloan filters i jjester et allboOSh . The X-ray data 
shows the unfolded PCA-l-HEXTE power-law model. 

The optical flux of the source is a significant frac- 
tion of its X-ray power, with monochromatic flux ratios of 
ALa(5500 A)/XLx(2 keV)=0.50. In terms of integrated fluxes, 

AI/a(5500 A)/L2-10 keV=0.46, or L4OO-7OO nm/i2-10 koV=0.30 

and Z/400-700 nm/ii-100 kev=0.07, respectively. Dereddening in- 
troduces systematic errors of a factor of «1.6 around 5500 A 

(§[I3- 



Both the r' and g' data show a mean flux that is slowly variable 
on timescales of several seconds or longer. Superposed on this are 
rapid flaring events. All this variability is statistically very signif- 
icant, as can be easily deduced by eye from the high labelled net 
count rates, as well as the flat comparison star light curve, shown 
for the r' filter. This comparison light curve is shown normalised 
to the median flux of GX 339-4 (and offset) for clarity; with this 
scaling, atmospheric scintillation equally affecting both the target 
and the standard star should cause variations of identical amplitude 
in both. It is seen from the comparison star light curve that the ef- 
fect of the atmosphere is minute, relative to the changes seen in the 
target, which must be intrinsic. 



4 RESULTS: TIMING 

4.1 Light curves 

Short 60 s sections of the net barycentered light curves from the 
best night (Night 1) are shown in Fig. [5] These are simultaneous 
for the three ULTRACAM filters, and the full band RXTE PCA 
GoodXenon events. The X-ray time resolution for this display is 
matched to the ~50 ms binning of the r' and g' filters, while the 
50 co-adds for each u' data point result in a time bin of «2.5 s. 



The novelty of our work is the fast time (<; 1 s) variability of 
the source at multiple wavelengths. Variability on slow er timescales 
has already been investigated in the literature (e.g. | Motch et al.l 
1 19821 ; ISteiman-Cameron et al.lll990l : iHoman et alj2005h . but iiToT 
der to ease comparison of results, we present in Fig.[6]l a longer 500 
s section of the simultaneous X-ray and r' light curves binned to 1 
s resolution (the variations in the other optical filters are similar). 
This section was chosen randomly with the only criterion that it 
contain several prominent X-ray flares, the strongest of which are 
marked according to the superposed flare decryption that we will 
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Figure 4. Broad-band SED of GX 339^, with simultaneous intrinsic op- 
tical (orange) and unfolded X-ray model (blue) spectral fluxes from Night 
1 . The green optical spectra represent the systematic uncertainty of the red- 
dening connection, computed for AAv=±0.5. 



next describe. The figure shows that there is no significant evolu- 
tion of the general light curve characteristics on this timescale. 

What is the maximum flare strength associated with the fast 
optical variations ? This can be answered by means of a su perposed 
shots analysis (cf. lNegoro et alJl200lLlMalzac et ai]|2003l Paper I). 
We selected flares with a peak strength of / times above the local 
mean in a running tm=32 s long section, and also required the flares 
to be maxima within a contiguous segment of ±ip=4 s. This assures 
that only distinct, significant flares are selected. The sections were 
all continuum-normalised, peak-aligned and averaged. The results 
for flares with /^2 are shown in the left panel of Fig.|7] The av- 
eraged flares have a peak value of close to a factor of 2 above the 
local mean, which means that strongest flares in both r' and g' are 
associated with flux increases of this magnitude. 

The plot shows the flares to be very narrow with widths of at 
most a few time bins. In order to quantify the typical time scales as- 
sociated with the optical flaring, we selected more common lower- 
amplitude flares with /^1.2 to which we fitted a simple exponen- 
tial time decay model of form Ae~*''^ . This is shown in the right 
hand panel of Fig.l?] where the profile of the average from flare su- 
perposition in both filters is plotted after removal of the normalised 
continuum of value 1. Overlaid are decay models fitted at times 
of <1 s from the average flare peak. We find r=0.11±0.01 s and 
0.09±0.01 s in r' and g' , respectively, implying that the fastest rise 
and fall timescales of optical flaring that we probe are only ~100 
ms. 

The variations in the two filters match each other well, and the 
peak to peak variability across the entire r' and g' lightcurves is 
a factor of 4. This can be discerned, for instance, in Fig.jS] where 
we plot every simultaneous flux measurement from the light curve 
of Night 1 in the two filters. The range of variations on both axes 
encompass a factor of 4. These lightcurves, when binned to the 
slower u data, also result in a good match to the u variations. 

Table[3]shows the average count rates in different bands. Also 
listed are the fractional rms variability amplitudes. These are mea- 
sured from the entire lightcurves and deno te the excess variati on 
above that due to Poisson uncertainties (cf. IVaughan et al.ll2003h : 
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Table 3. Light curve properties. R represents the net counts per second 
(rounded off to the nearest 10 ct s^^ in the optical), hi X-rays, this is stated 
per active PCU in the PCA instrument, rms is the fractional variability am- 
plitude above Poisson uncertainties (Eq.[T). Typical \u error on this value 
is <0.02. Light curves with time resolutions matching the best optical time 
resolution on each night were used for these measurements in both opti- 
cal and X-rays; only the u' observations are 50 times slower (see Table[T). 
Note that non-photometric weather causes the mean optical flux to decrease 
significantly on Nights 2 and 3. 




200 300 

Time (s) 

Figure 6. Simultaneous 500 s section of extracted net barycentered 
lightcurves from Night 1 in ULTRACAM r' and X-ray PCA full band data, 
binned to 1 s time resolution. X-rays flares with /^2.4 are marked. 



where I denotes the mean count rate of the light curve I in any band, 
var is its raw variance, cTorr the individual errors on the count rate 
measurements and cr|rr the mean square error. 

The X-ray light curves were extracted with a time bin of 2~* s 
from the full PCA band. In addition, light curves in two other bands 
were extracted for comparison: ~2-5 keV (PCA channels 5-11) 
and ?5i5-20 keV (PCA channels 12-47). The X-ray count rate is 
much lower, and fractional variability higher, as compared to the 
optical. The strongest optical variability occurs in the r' band, at 
about 15 per cent rms across the full light curve. This red rms 
variability is perceptibly stronger than in g , as can also be seen 
from the slope of the fit to the distribution of instantaneous fluxes 
in Fig.IH 



4.2 Power spectra 

The power spectral densities (PSDs) are computed as the modulus- 
square of the discrete Fourier transforms of the light curves as 
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Figure 5. Simultaneous section of extracted net baiycentered liglitcurves from Niglit 1. Tlie ULTRACAM r' and g' data liave a time resolution of 50 ms, and 
the X-ray PCA full band data (summed over two active PCUs on this night) is matched to this for comparison. The u' data is 50 times slower The lower dashed 
curves in all optical filters are for the comparison star and normalised to the mean flux of GX 339^ in each filter (with an offset downwards for clarity), so 
any atmospheric scintillation will appear as equal amplitude modulations for both target and comparison. 
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Figure 7. Peak-aligned optical flares properties, normalised to the local mean flux level, (heft) The strongest flares, selected with /^2. The lower (green) 
histogram is the average for g' and the top (red) histogram is for r' (offset by -hO.4 for clarity). In both cases, the peak has a value of about twice the local 
mean level. (Right) The superposition of more typical flares selected with /^1.2 is shown, overplotted with best-fit exponential decay profiles with a decay 
time of 0. 1 1 and 0.09 s in r' (red) and g' (green) respectively. 
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Figure 8. Instantaneous light curve count rate measurements for all 68317 
time bins from 56 mins of data in r' vs. g' on Night 1, normalized to the 
mean flux in each filter. The 1 : 1 line is the dashed light blue one, and the red 
continuous line is the best-fit to the linear model ofy = mx + c obtained 
by the ordinary least square bisector method of llsobe et al J ll990), with 
m=0.871±0.002 and c=0.129±0.002. The gradient is flatter than the 1:1 
line, which means that flaring is stronger in r ' . 



with the frequencies themselves being discretised at N/2 evenly- 
spaced values from l/NAT to the Nyquist fre quency 1/2AT, 
where AT is the time bin (cf. lVaughan et alj2003h . 

The X-ray PSDs were computed from net hghtcurves with 
time resolution of 2~* s over the full PCA energy range. The data 
from each night was split into segments of 256 s, and the power 
spectra were computed using standard XRONOS procedures, and 
cross-checked with a custom routine coding Eq. |2] and also used 
for the optical PSDs below. In this case, I refers to the mean count 
rate in any given segment under consideration. The PSDs were av- 
eraged across the segments and rebinned in logarithmic frequency 
bins by a factor of 1.1 in order to improve signahnoise, especially 
at high Fourier frequencies. The optical PSDs were computed using 
an identical procedure, on the fastest time resolutions available for 
each filter and each night. A rebinning factor of 1.05 was sufficient 
for all data on Nights 1 and 2. On Night 3, no acceptable solution 
was found for the PSD fits that we will describe shortly, likely as 
a result of bad weather systematics. For this night, we rebinned by 
a factor of 1.1. The constant A in Eq.|2]was determin ed following 
usual rms-squared standardisation A=2AT/P N (e.g. Ivan der Klisl 
ll997l ; lBelloni&Hasinge3ll99d : IVaughanetal.ll2003l) . which nor- 
malises the PSD to the fractional light curve variance per unit fre- 
quency interval (above Poisson uncertainties if white noise is sub- 
tracted, as described below). The PSDs are shown in Fig.|9]for each 
night, in both optical and X-rays. 



This proved to be an excellent measure of the constant noise level in 
the X-ray PSDs, as is readily apparent in Fig.|9] It also showed PCA 
dead time corrections to be unimportant in determining this noise 
level for th e low count ra t e relev ant here; this was confirmed by us- 
ing Eq. 4 of lNowak et al.l ( ll999ah to determine the correction factor. 
Poisson uncertainties on the gross (source-l-background) count rates 
are the only important contribution to the noise in X-rays. In the 
optical, additional contributions to the variance are CCD readout 
noise and dead time. The read noise contribution is automatically 
accounted for in the photometric errors returned by the pipeline. 
The presence of dead time means that any fast flux variations that 
occur within the dead time interval are not sampled and cannot be 
distinguished from source variations on (longer) sampled times. Ef- 
fectively, some power is 'aliased' down to fre quencies below th e 
Nyquist sampling frequency of l/2rcycie (e.g. Ivan der Klisll989l) . 
The high duty cycle of our observations (Table 1) means that alias- 
ing will be small and can be approximated by incre asing the noise 
level by the ratio of Tcydc/Tcxp; see Eqn. (A3) of IVaughan et al.l 
2003. The resultant corrections increase the noise power by « 2-3 
per cent between the nights (cf. § 12. U and have been applied. Fi- 
nally, photometric errors from the comparison star against which 
relative photometry is performed are included in the final noise 
level via Gaussian propagation of the errors to crlrr in Eq. [3] The 
noise levels are indicated on Fig.|9|and show all optical data to lie 
above the white noise level up to the Nyquist frequency. Several 
tests of the photometry and errors were carried out to confirm this, 
and the effect of any possible systematic misestimate in noise levels 
is discussed alongside other results presented below and in the fol- 
lowing sections. In any case, the optical noise levels are generally 
low because of the high optical count rates (TablejS); consequently, 
their effect on the mean PSD levels is much smaller than in X-rays 
over the full frequency range probed. The increasing optical noise 
level on Nights 2 and 3 can be explained by a significant decrease 
in detected source count rates through high and strongly variable 
transparency and seeing conditions. This also degraded the quality 
of the slower u data, for which only the Night 1 data are presented. 
The results after noise subtraction are displayed in the bottom 
row of Fig. |9| for each night of data separately. The X-ray vari- 
ability power is larger than that in the optical on all nights. Noise 
dominates the X-ray PSDs at Fourier frequencies of >10 Hz, with 
the result that using light curves with longer binning (e.g. 50 ms 
matched to the optical) produce identical PSDs to those made from 
the above finer time binning. The optical power spectra also show 
band-limited noise continua up to several Hz at least. Above ~0.1 
Hz, they exhibit strong curvature to high Fourier frequencies. The 
limiting Nyquist frequency decreases from Night 1 to Night 3 as a 
result of the increasing time bins with worsening weather. Around 
~0.05 Hz, there is a local excess on all nights, whereas at the low- 
est frequencies, the continuum slope seems to exhibit a change. 
Despite the changeable weather, the overall PSD structures on all 
three nights match. 



4.2.1 Noise levels 

Poisson uncertainties introduce a constant white noise level, which 
becomes increasingly important at high Fourier frequency. For in- 
dividual (mean) count rate errors (Jcrr (o^), the noise for rms- 
squared normalised PSDs given by 



4.2.2 Model fits 

In order to parametrise the PSD shapes, we fit them with Lorentzian 
profiles. Given the low source X-ray count rate, the relatively short 
total observing times, and the fact that the PSDs appeared simi- 
lar on all nights, we averaged all three nights in order to charac- 
terise the functional form. This was done separately for the full 
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PCA band, as well as the 2-5 keV and 5-20 keV bands. Three 
Lorentzians were found to produce a good fit in each case, as shown 
in Fig. [To] These are all broad (or zero-centred) Lorentzians de- 
scribing the broad peaks at characteristic frequencies (i/max) of 
around 0.02, 0.1 and 2 Hz, consistent across the three energy bands. 
The measurements are listed in Table [4] No quasi-periodic oscilla- 
tion (QPO; defined as having a quality factor of Q > 2) is found, 
with the narrowest Lorentzian having Q~0.7. The Lorentzian pa- 
rameters - r, t^o. A, z^max and Q ~ are defined identically to 
lBeUonietal]J2002h . 

In the optical, the PSDs were similarly fit with multiple 
Lorentzians. All PSDs clearly contain a low-frequency QPO and 
the PSDs covering fast timescales (i.e. the r' and g' data) also de- 
mand extra curvature at high Fourier frequencies (Fig.|9j. This is 
especially clear for the Night 1 data which provide the best dataset 
in terms of probing the broadest range of timescales, and observed 
under the best weather conditions. A minimum of four Lorentzian 
components are required to describe this - one for the QPO and 
three for the broad-band continuum humps (Fig. II II Table|5]l. We 
also tried different models, including an ex ponential cut-off power- 
law, as well as the bending power-law of iMcHardv et alj ( 12004) . 
for the PSD fall-off above a Fourier frequency of ~1 Hz. These 
can roughly describe the curvature over restricted frequency ranges, 
but require further modification such as additional Lorentzians or 
multiple bending sub-components, for a reasonable fit. Given our 
relatively short dataset, we chose to retain the multiple Lorentzian 
component fits described above as a simple parametrisation of the 
PSD, and not necessarily as an accurate physical description of the 
underlying variable components, though as we will describe in the 
following sections, this parametrisation does provide interesting in- 
sight into the X-ray vs. optical coherence and time lags. 

Note that we have not included the lowest and highest Fourier 
frequency bins in the optical fits presented, because the presence of 
new components are indicated at both ends which cannot be prop- 
erly modelled. We have kept the fits to below 8 Hz. Including the 
bins up to 10 Hz changes x^ significantly; e.g. Ax^=+20 for an 
extra five degree of freedom for the r' fit. The presence of vari- 
ability noise above our Nyquist sampling frequency is important, 
and should be probed in future faster observations. A lower limit 
on the strength of the high frequency component may be obtained 
by extending the fits shown in Fig. \TT\ and measuring the residu- 
als, which gives an rms of at least 0.5 per cent, over only 8-10 
Hz. We also note that the fit is robust to small systematic changes 
in the white noise level. We tested this by subtracting off a higher 
noise obtained from fitting a constant (per unit Hz) to the PSD at 
the highest Fourier frequency bins; the relative normalisation of the 
individual components does vary, but not the overall fit. It is also 
worth noting that in the X-ray PSD (Fig.[TOl top panel), there is an 
apparent excess above the three Lorentzian fit at the high frequency 
end (centred at ~9.4 Hz). But fitting an extra Lorentzian to this re- 
sults in a Ax^ of only -2 for 3 extra degrees of freedom; i.e. this 
excess is not significant in the present data. 

We carried out similar multi-component Lorentzian fits to the 
r' and g' data for Nights 2 and 3, which cover a narrower range 
of Fourier frequencies. The PSDs for these two nights are pre- 
sented in Figs. 1121 and 1131 and fit parameters listed in Table [5] 
Three Lorentzian components were sufficient for the data on Night 
3 (which cover the narrowest frequency range, and have been re- 
binned by a factor of 1.1 as mentioned before). The characteris- 
tic frequencies of the QPO and the component covering the peak 
around 1 Hz show very mild variations across the nights, with an 
apparent increase on Nights 2 and 3. This could simply be a re- 



sult of the fact that we do not probe the high Fourier frequency 
Lorentzians well on these nights, which forces the fit peaks to shift 
slightly in this direction to compensate. In any case, there is con- 
sistency between these fit values at the 4 a level across the nights. 
At low Fourier frequencies below ~0.02 Hz, both Nights 2 and 3 
suggest the presence of some small excess power above the Night 1 
variability power levels. This may be indicative of incomplete pho- 
tometric correction when dividing by the comparison star, because 
any residual transparency variations are likely to be more important 
on long timescales (low Fourier frequencies). Additional system- 
atic residual effects are apparent on the worst Night 3, where the 
QPO appears to be broader (A~0.01 1 in r', compared to a value of 
0.004 on the other two nights), and is not well described by a single 
Lorentzian. 

Within the limits imposed by the above minor caveats, the 
PSDs may be considered to be stable across the nights, and for the 
physical interpretation presented in the remainder of the paper, we 
will mainly refer to the Lorentzian fits for the best Night 1 . A more 
detailed analysis of the PSDs will also be presented in future work. 



4.2.3 The QPO 

The QPO lies at a Fourier frequency of ~0.05 Hz (on all nights) 
and has a Q value of 5-6 on the best Nights 1 and 2; Table[5](. The 
u' data covers only low Fourier frequencies, but a simple fit to this 
clearly shows the QPO to have similar properties in all three fil- 
ters (Fig. II lb. The QPO is not an instrumental artifact, otherwise 
it would have introduced oscillations in the comparison star light 
curve as well. It is clear from Fig. [5] that this is not the case. We 
also confirmed this by extracting optical PSDs for comparison star. 
Additionally, there is no known mechanical or electronic cycle ei- 
ther in the instrument, or in the telescope which occurs at these 
frequencies. 

We also checked the detectability of any low frequency QPO 
in X-rays by adding a model Lorentzian to the PSD fit of Fig. [TO] 
with a fixed i^max and Q identical to that measured in the optical, 
and obtained a 90% (Ax^=2.71) upper limit on its integrated frac- 
tional rms of ''x-ray ~0.075. This is about 2.5 times larger than 
^Optical ~0-03 (see Table[5](, so we cannot rule out an X-ray QPO 
identical to that seen in the optical. 



4.2.4 Previous PSD measurements 

How do these PSDs compare with previous measurements? In X- 
rays, there have been several studies. For example , our Lorentzian 
fit is similar to that found by iNowak et al.l j 19991) for their obser- 
vation 5, whic h had the faintest source fluxes (_F3_9=2.5x 10~^" 
erg s^^ cm^^ I Wilms et al.lll999r . only about twice as bright as in 
our observation). In particular, the characteristic frequencies of the 
two strong broad band features match closely. In the optical, previ- 
ous measu rements hav e been published below a Fourier frequency 
of 1 Hz bv 'Motch et a l. ( 1982, 1983) durin g a 1981 low/hard state 
(bright in both optical and X-rays), and by Motch et al.l ( 119851) dur- 
ing a 1982 X-ray off state. The hard state observation showed a 0.05 
Hz (20 s) QPO feature, very similar to our observation. The QPO 
was found at a higher frequency (0.14 Hz, or 7 s) in the off state. As 
these data were obtained with a completely different observational 
setup, this also rules out any instrumental artifice for our QPO de- 
tection. The high frequency power decreased rapidly in all these 
observations (e.g. with a slope of -1.6 in 1981), whereas we find a 
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gentler initial fall-off; fitti ng the 0. 1-1 Hz PSD with a single power- 
law for comparison with iMotch et alj J1982r) gives a slope of ~- 
0.5. We note that the bright hard state observation o f 1981 showed 
a strong QPO in X-rays as well JMotch et al.lll985r) . which is ab- 
sent in our (fain ter) case. Other studie s incl ude very fast white light 
observ ations bv llmamura et alj ( 1 199(1) and lSteiman-Cameron et al.l 
( Il997r) . though with no simultaneous X-ray data. The last authors 
detected strong continuum curvature around 1 Hz and a 0.06 Hz (16 
s) QPO during 1996, suggesting some reasonable agreement with 
our r' and g' PSDs. 



4.3 Auto-correlations 

The auto-correlation function (ACF) represents a convolution of a 
light curve with itself, and gives a measure of the effective coher- 
ence timescales in the data. For a regularly-sampled time series /, 
a simple (but biased) estimator of the ACF at any given lag r, may 
be computed as 



A^ir^^^^^^^^-"^^^ 



I] 



Y.l,W. 



lY 



(4) 



where the lags are discretised so that tj^ =ti + tj . If Poisson noise 
is strong, auto-correlation of the errors causes a prominent, narrow 
spike to appear at zero lag (normalised to 1), thereby suppressing 
the relative contribution at all other lags. Random errors are un- 
correlated for all non-zero lags. A simple way to correct for this 
component is by use of the Wiener-Khinchin theorem, which states 
that the power- spectral density of a random process and its auto- 
correlation function are Fourier pairs. This allows removal of the 
Poisson contribution by correcting the PSD for white noise power, 
followed by an inverse transform. 



ACF(t) 



E 



PSD{f) - n]e 



27TifT 



(5) 



where the symbols have been described in the previous section, and 
the ACF is defined to be normalised to 1 at zero lag. Removal of 
the noise (n) is crucial for the X-ray data as expected from the 
low count rate. The resultant noise-corrected full band PCA ACF 
lies above the uncorrected one by a factor of ~5.8 (at non-zero 
lags) on Night 1 , for example. This low rate limits the detectability 
of real, high frequency variable components, and we will consider 
the implications of the presence of any such fast variability in the 
Discussion section. On the other hand, white noise is only a minor 
correction to the fast optical PSDs, as discussed in the previous 
section. This is a result of the high count rate, which means that 
n lies much below the mean source variability power over the full 
Fourier frequency range (see Fig. [9] and Table [3]l. The maximum 
correction necessary is for the slower u' data where the count rate is 
lower, but even in this case, the PSD lies above the noise level at all 
frequencies so n does not dominate. The resultant noise-corrected 
ACFs in this case lie above the uncorrected ones by only a factor 
of 1.4. The corresponding correction factor is much smaller in the 
r' and g' bands, at 1.1 and 1.2 respectively. 

The noise-corrected ACF results for the three nights of data in 
the different bands are shown in Fig. [14] For exact comparison, the 
X-ray ACFs are constructed from light curves binned to the fastest 
optical time resolution on each night. On all nights, the fast op- 
tical data (r' and g') have a narrower central peak relative to the 
X-ray ACF. The u' light curve is sparser, and the effect of white 
noise comparatively larger on its PSD. Yet, the u' data show simi- 
lar characteristics, with a core marginally narrower than the X-ray 



ACF computed from an identically binned X-ray light curve with 
2.5 s time resolution. Such 'slow' ACFs for all bands are shown in 
the bottom panel of Fig. [14] 

Another feature in all the optical ACFs is the positive hump 
around lags of 20 s, a result of the optical QPO discussed in the 
previous section. Some weak humps are also present in the X-ray 
data between lags of 15-25 s, though their significance is small. 
These can be quantified simply by first computing the standard de- 
viation of ACF measurements from independent light curve seg- 
ments, and using this standard deviation as the error estimate at 
each time lag bin. Then, the positive ACF signal over the 15-25 
s is coadded and compared to the propagated error on this coad- 
ded value. This results in significance estimates of 3.8ct, 2.5(j and 
<la for the X-ray humps on Nights 1, 2 and 3 respectively. Some 
of these may be related to the broad Lorentzian 2 component with 
J'max^O. 1 Hz found in the X-ray PSD fit (Table |4]l; as this is not 
a QPO, no strong ACF excess is expected. Carrying out the same 
computation in the optical shows all the ~20 s lag humps to be 
highly significant (>>10(t). 

There is some slight apparent evolution of the ACFs between 
the various nights. For instance, the optical ACFs on Night 3 are 
marginally broader, and the X-ray ACF narrower, than in the other 
nights. Differences in optical weather conditions, the minimum 
time resolutions, and the overall short observation lengths, must be 
kept in mind (Table[T]l, so detailed inter-night comparisons are more 
involved, and not considered here. Regardless of these changes, 
the main ACF structure within lags of ±10 s possess several com- 
ponents on all the nights, with a very narrow central core and a 
broader component on longer lags. For Night 1, the core has an 
exponential decay time of close to 0. 1 s, a manifestation of the su- 
perposed flare analysis shown in § 14.11 In X-rays, this fast decay 
time is closer to 0.2 s. Lastly, we also note that small systematic 
changes in the optical PSD white noise level do not affect the over- 
all comparisons presented herein, with the fast optical ACFs still 
being much narrower than in X-rays if a higher optical PSD noise 
component (as detailed in the previous section) is subtracted. 



4.4 Optical vs. X-ray cross-spectral behaviour 

4.4.1 Time-domain cross-correlation functions 

JMotch et al.l ( Il983l) uncovered a weak anti-correlation between the 
optical and X-ray flux around the time of a source transition from 
the low to the high state. In our simultaneous VLTAJLTRACAM 
and RXTE observations, we confirmed an r' anti-correlation on 
timescales of a few seconds, both leading and lagging the X-ray 
peak by a few seconds, and also reported the presence of a sharp 
positive correlation spike lagging X-rays by only fractions of a sec- 
ond (Paper I). A superposed shot analysis confirmed the broad op- 
tical anti-correlation and also revealed fast optical spikes lagging 
the X-rays by 150 ms. 

In Fig. [15] we now present the CCFs for all three filters, on 
the various nights. The CCF represents the convolution of one light 
curve with respect to the other. For our discrete time series, we used 
a standard IDL coding for the time-domain convolution, which is 
estimated as 



CCF(t,) = <^ 



^"^^M^(tj.)-^][°(t.)-°] 

^"~^'b(i.)-S^l[°(t,,)-°] 



r < 



T ^ 



(6) 
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Figure 9. Power spectra of GX 339^ for Nights 1-3 (from Left to Right). The X-ray full band PSDs are in black, while the optical ones are in red, green 
and blue for r', g' and u', respectively. The top row shows the PSDs for the net light curves without subtraction of the white noise contribution, which is 
indicated by the dotted horizontal lines in the respective colours. Power on the y-axis is per unit Hz, with the standard rms-squaied normalisation. In these 
units, the white noise is a constant, and the integral of the PSD over frequency results in the square of the fractional rms amplitude. The bottom row shows the 
noise-subtracted PSDs, now in units of frequency X Power. 



where tj^ =ti + Tj. x and o represent the first and second light curve, 
in this case the X-ray and optical ones, respectively. The lag Tj 
refers to that of the second light curve, with respect to the first, 
and is itself discretised. Light curves which are related only by a 
linear transfer function (e.g. a simple time delay or a multiplicative 
scaling factor) would give a peak CCF value of 1 at the relevant lag 
time. 

The above formula may be used if the two light curves have 
been extracted on an identical time baseline. This was possible be- 
cause of the event mode of the RXTEFCA, with every single event 
being recorded. X-ray light curves simultaneous to the optical were 
obtained by binning the fast (2^* s) barycentred X-ray light curves 
on the optical baseline, separately for each filter and night. We also 
checked that direct FTOOLS extraction of X-ray light curves with 
time bins matching the optical ones, followed by minimal interpo- 
lation to the common time baseline, produced unchanged results. 
Such a methodology is standard procedure in cross-correlation 
AGN light curves (e.g. lOaskell & Petersonlll987h . and has been 
shown to produce reliabl e results in the limit of well-sampled data 
d White & PetersonI 19941) . As an extra check, we also computed the 
discre te correlation function (DCF) described bv Edelson & KrolikI 
( 119881 their Eq. 2). The DCF does not require any interpolation 
of the data itself, and may be computed using datasets with differ- 
ing time resolutions. Rather, a correlation is computed between all 
pairs of data in the two light curves, and a time lag is attached to 
each pair. The final DCF is then measured by averaging all the pairs 
that contribute to a binned time lag series. As expected, the match 
in the CCF shapes between the two methods was excellent for a 
variety of time lags and X-ray light curve time resolutions. Note 
that no additional correction for measurement error was applied 
(cf. IWhite & Petersonl ll994f) : in any case, this would only affect 
the DCF normalisation. 

Each simultaneous light curve was then split into segments 



256 s long, which were cross-correlated and then averaged to ob- 
tain the final result. The uncertainties on the average at each lag 
time bin are estimated by error propagation, using the standard de- 
viation of the CCF among these segments. The main CCF structure 
within absolute lags of several seconds is reproduced on all the 
nights, with the r' and g' results matching closely in both strength 
and shape. The most significant other feature is a positive correla- 
tion hump around lag=+10 s. Using the coadded signal methodol- 
ogy used for estimating hump significance in the ACFs (i) l4.3t . this 
time between optical lags of +5 and -1-15 s, results in a signahnoise 
of 10-30 for this feature between the various nights and the r' and 
g' filters. Further structures at lags greater than -1-20 s and below -10 
s are also present, though these change dramatically in strength and 
width between the nights these are not considered to be sigificant 
in the full data. The main features have been noted in Paper I with 
regard to the r' CCF and the overlaid and zoomed-in CCFs from all 
nights presented in Fig. 1 of that paper make immediately apparent 
the significance of the main features being discussed (note: in that 
work, the length of the time segments was not fixed to 256 s, which 
explains minor differences on long lags with respect to Fig.|15|here. 
but all the main features appear in both cases). Again, stringent con- 
straints on the evolution of the CCF between the nights must await 
better datasets, given the changing weather and time resolution of 
the present observations. In particular, the prominent variation in 
the CCF on Night 3 on times of several tens of seconds is an arti- 
fact of highly variable transparency. 

The u' data are much slower; at a time resolution of 2.5 s, each 
256 s-long u' band segment has only 103 time bins (as opposed to 
5120 bins for the other filters). In Fig.[T5] we compare the X-ray 
vs. It' CCF from Night 1 with slow X-ray vs. r'/g' CCFs, with the 
latter computed from light curves heavily binned by a factor of 50. 
Two features stand out from this figure:- i) The u' shows an anti- 
correlation at small negative lags (~-5 s) as in the other filters. The 
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Figure 10. Lorentzian fits to the X-ray PSDs in the full PCA (Top), the 2-5 keV (Bottom left) and the 5-20 keV (Bottom right) bands. The X-ray data are from 
the average of all three nights. See Table|4]and text for fit parameters and details. 



formal significance of this feature is only 2cr, but this is a result 
of the large reduction in light curve time bins (and hence indepen- 
dent segments) as compared to the fast data for r' and g' . The fact 
that the location of this trough agrees in all filters means that it is 
real, ii) The prominence of the ~-f10 s lag strength increases in u 
(reaching a maximum value close to 0. 15 over the plotted lag range) 
accompanied by a decrease in the fast sub-second lag features ap- 
parent in the other filters. The difference between the CCFs com- 
puted from fast and from slow light curves may be understood as a 
result of the very different widths of the respective CCF structures - 
the fast and strong peak is confined to very small lags, whereas the 
~10 s structure is broad and is therefore enhanced when binned. 

Despite the slow u time resolution, it is possible to check 
for the presence of any fast variations correlated with X-rays. We 
did this by using the slow u' light curves, and the fast X-ray light 
curves sampled at 50 ms resolution, and (because of the differing 
time resolutions) computing the DCF between these two. The result 
is shown in Fig. [161 where the final DCF is also shown with bins 
corresponding to the X-ray time sampling. The errors are computed 
according to the prescription of lEdelson & KrolikI jl988r) , with the 
uncertainty on the DCF depending simply on the scatter of the in- 
dependent lag values that contribute to each bin. Other than the 
dominant slow component increasing toward 10 s (cf. last panel of 
Fig.|15t, there is a weak excess on sub-second lags, just as for the 
r' and g' filters. This feature appears to be broad (extending over 
the central ~1 s; note that we cannot claim a significant detection 



of a 150 ms lagged peak itself), but this is not unexpected given 
the large u time bins. The overall significance of the broad fea- 
ture can be assessed by co-adding bins to obtain a mean DCF and 
comparing this to the propagated error on the mean (as for the ACF 
humps), from which we find a significance of 5.8(7 over the lag 
range of -1 to 1 s. Alternatively, fitting a simple line profile (e.g. 
a Gaussian) and comparing this to a constant model fit returns an 
improvement in the former case at a level of greater than 99.99 per 
cent, according to an f-test. Similar excesses are also detected in 
DCF computations over independent light curve segments. Thus, 
the excess is real. This exercise shows that the overall fast vari- 
ability patterns appear to be present in all three filters (though with 
different strengths). 

What about the behaviour of one filter vs. another? As is seen 
from good match between the various optical light curves, the CCFs 
between the different filters do not reveal any strong correlation 
characteristics. The r' vs. g' CCF is shown in Fig. 1171 and is seen 
to be largely symmetric with no lag of the peak. Intriguingly, there 
is a very small asymmetry discernible in the sense that the CCF is 
slightly stronger on positive r' lags. This effect is seen on multiple 
nights, and so is probably real, but minute. We hope to observe the 
source with finer time resolution in the future, and defer discussion 
of this for now. 

We also investigated the X-ray-only behaviour. We used the 
RXTE light curves extracted in two separate bands: 2-5 keV and 
5-20 keV (see Table [3} and cross-correlated them in an identical 
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Figure 11. Optical PSD fits to Night 1 data. Four component model fits to the r' (Top Left) and 'g' data (Top Right). The bottom panel is a three Lorentzian fit 
to the u' . See Table|5]and text for fit parameters and details. 
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Figure 12. Four component model fits to the Optical PSD fits to Night 2 data, for r' (Left) and ' g' data (Right). See Table[5]for fit parameters. 



manner as described above. The results are also shown in Fig. 1151 
with a positive lag referring to the hard energies lagging soft pho- 
tons. In contrast to the asymmetric X-ray vs. optical CCFs, the soft 
vs. hard X-ray CCFs are largely symmetric and show no obvious 
peak lag. This can also be seen in Fig.ll7lwhere the CCF computed 
on the fastest X-ray time bin of 2^* s is shown zoomed-in. There 



are no obviously significant structures to within the sensitivity lim- 
its. 



4.4.2 Coherence and time lags 

The time domain cross-correlation functions can be de- 
composed into their Fourier components - coherence and 
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Figure 13. Four component model fits to the Optical PSD fits to Night 3 data, for r' (Left) and ' g' data (Right). See Table[5]for fit parameters. 



phase lags j jenkins & Watts' '1969"; 'Bendat & Piersoll 1 19861 : 
IVaughan & No wak 1997; Nowaketal. 1999a) . In order to com- 

f ute th ese, we follow the recipe described by IVaughan & NowakI 
19971), coded with a combination of custom Fast Fourier Trans- 
form routines and a modified version of the publicly available 
package IDL.EXTRAC'o The simultaneous lightcurves were split 
into many segments as with the PSD measurement, and the cross- 
spectrum between each X-ray and optical segment was computed. 
A coherence function is then computed from the average of many 
segments as 



lU) 



\{C)\' 



{\S.?){\So?) 



(7) 



where x and o represent the X-ray and optical light curves, X 
and O their Fourier transforms and C the complex-valued cross- 
spectrum between these two C = X*0 sX any Fourier frequency 
/. Poisson noise correction is included via the S and rixo terms. 
jSp represents the PSD at frequency / after correction for 



white-noise: i.e. \Sx\ ~ \X\ 



liVJ"and|Sor = |0|- 



\No\ 



with TV being the Fourier transform of the noise component. The 
PSDs are typically written in Leahy normalised units so that the 
ensemble noise level is 2. Uxo is a combination of signal and noise 
power terms: nL = {\Sx\''\No\' + |iV.|'|5op + \Nx\''\No\^)/m 
for m independent light curve segments and frequency bins which 
are averaged over. For all frequency bins satisfying \S\^ /\N\^ > 
cr/-^/m (for both light curves), |{S^So)|^ > crn^o/-ym and 7^ > 
o"'^xo/(|-''^P|0|^)i the final corrected coherence can be analytically 
estimated from Eg. [7] above. The u ncertainty is computed as in 
Eqn. 8 of IVaughan & NowakI <1997h . This effectively refers to the 
Gaussian limit and means that it is only valid for significant coher- 
ence and PSD power levels, for which we used a threshold value of 
(T=3. The computation was performed on geometric ally-rebinned 
(in frequency) cross- and power spectra. 

The phase lag ((j>) is simply the phase angle of the complex- 
valued cross spectrum defined over the range — tt to tt. If one takes 
7?aw to be the raw coherence function before noise correction so 
that7r\^,, = |{C)p/({|X|^){10p)), then the phase uncertainty is 
given by 54>{f) ~ v (1 ^ 7law)/27?aw"i. We refer the reader to 



the works cited at the beginning of this section for full details on 
the theory. 

* http://idlastro.gsfc.nasa.gov/ftp/contrib/rxte/ 
© 2010 RAS, MNRAS 0Q0.fTH27l 



The results are shown in Fig. [Ts] As the coherence function 
requires averaging over many segments, the data from the best 
two nights (1 and 2) have been averaged for the high time reso- 
lution r' and g' light curves over the common frequency range, and 
heavily rebinned in frequency. The data from Night 3 (not shown) 
have a similar pattern, but with a lower mean coherence and larger 
(weather-related) scatter. There are several small differences in de- 
tail in the PSD fits on the two nights (^ I4.2.2l and Figs.[TT][l2ll, but 
given the wide Fourier frequency bins that we use, the small PSD 
changes between the nights are unimportant. 

At high Fourier frequencies, the coherence function shows a 
broad profile with a peak around 0.7 Hz, and a maximum value of 
~0.1 reflecting the comparatively- weak CCF strength. Its value de- 
creases at both low and high frequencies, but remains significantly 
above zero over at least the range of 0. 1-5 Hz. There is good overall 
agreement between the two filters. The lower Fourier frequencies 
exhibit a rise again, and reach a maximum coherence of «0.4 in 
the blue u' filter. Around 0.1 Hz, the coherence values are approx- 
imately consistent between the filters. 

The argument of the complex-valued cross-spectrum is the 
phase lag between the bands at any given Fourier frequency. These 
lags are shown in the third panel of Fig.[T8] One clear pattern is the 
smooth increase from ~0.2 Hz up to ~2 Hz. The phase lag is only 
defined between -n and n, so once it crosses the upper limit, the 
lag jumps down to -n. This holds true for all 27r multiple intervals 
of the phase. Thereafter, it should continue to increase again if the 
trend continues to higher Fourier frequency. The lags rising from 
— TT above ~3 Hz indicate that exactly such behaviour is occurring. 
The corresponding time lags at frequency / can be obtained by di- 
viding the phase lag by 2nf (Fig.[T8] Bottom). The rising phase 
lags correspond to a flat time lag distribution, with value of ~150 
ms at the maximum of the coherence function, exactly the peak 
time delay of the fast CCFs (Fig. 1151 Paper I). The behaviour of 
the coherence and lags at low Fourier frequencies is complex, with 
negative lags over ~0.08-0.3 Hz, a phase lag of ~7r/2 going down 
to 0.01 Hz, and another reversal with phase lag close to zero in 
the lowest bin. All these are associated with longer time lags over 
~l-20 s (mainly scattered around 10 s). 

The X-ray only behaviour of the coherence and lags (not 
shown) was found to be typical for XRBs. Computing the co- 
herence and lags between the PCA soft (2-5 keV) and hard (5- 
20 keV) light curves gave results very similar to the GX 339-4 
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Figure 14. The auto-correlation functions for Niglits 1-3 from top, going 
down. In each case, the black histogram is that for the full band X-rays 
(summed over all the PCUs available on each night), red for r' and green 
for g'. The effect of white noise has been removed in each case (see text). 
The last panel shows the ACFs for the photometric Night 1 , but computed 
from light curves matching the u' time bins of Ri2.5 s. In this case, the X- 
ray ACF is shown as the black dashed histogram, and the u' one is the blue 
histogram. 
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Figure 15. The cross-correlation functions (CCFs) between various bands, 
on Nights 1-3 (from top, down). The red and green curves show the X- 
ray full band PCA vs. the r' and g' CCFs, respectively. Positive lags mean 
that the optical is delayed with respect to X-rays. The black curves show 
the X-ray 2-5 keV vs. 5-20 keV CCFs (at the same time resolution as the 
optical, and normalised to 0. 1 for comparison), with delays referring to the 
5-20 keV band. For a zoom-in on the fast peak delay at 150 ms, see Fig. 1 
of Paper I. The bottom panel shows slow X-ray vs. optical CCFs, computed 
from light curves binned at 2.48 s, to match the u' data; in this case, the blue 
histogram refers to the X-ray vs. u' CCF. The error bars come from eiTor 
propagation using the standard deviation among independent light cu rve 
segments. © 2010 RAS, MNRAS 000,[Tlt22] 
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Figure 16. Zoom-in of the discrete conelation function (DCF) computed 
from fast X-ray data (sampled at 50 ms) vs. the slow u' data (2.5 s resolu- 
tion). The lag refers to the u' . A positive correlation around zero lag is seen, 
in addition to the rising, slower ~10 s component of Fig. 1151 The r' and g' 
fast CCF peak lag time of 150 ms is denoted by the arrow. 



study of iNowaket alj ( Il999l) . but now extended to lower source 
fluxes by about a factor of 2 - we found a high coherence value 
of ~0.9-l over much of the Fourier frequency range, and a log- 
arithmic time lag decline (within the errors) to high frequen cies 
(e.g. . Miyamoto et al. 1992; Vau ghan & Nowak.l997, ; ,Nowak et^ 
ll999bl : lKotovetalj200ll) . 
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5 DISCUSSION 

The X-ray spectral properties of GX 339-4 during our observa- 
tion are typical of th o se of the source in its X-ray low/hard state 
( lTomsicketaI]|2008l : I Wilms et all 1 19991), as are the X-ray PSD 
shape and the hard vs. soft band time lags I INowak et al.ll999n . The 
X-ray rms v ariability amplitude (Table [3) is high, but not unprece- 
dented (e.g. iBelloni et alj|2002!) . The optical flu x of the source is 
also w ithin the range observed during this state JBuxton & Bailvnl 
l2007h . The multi-filter optical PSD analysis and fits across three 
decades in frequency, and extending to ~10 Hz, is the most de- 
tailed so far for this source. The detection of an optical low- 
frequency QPO is not new (cf. J Motch et alJI 19831 ; llmamura et al.l 
I1990I : ISteiman-Cameron et al.lll99(J) . but is now confirmed in sev- 
eral filters. 

The origin of the source emission at multiple wavelengths 
remains uncertain, and the disc, relativistic jet plasma, or non- 
relativistic cyclosynchrotron from a cor ona are all proposed 
as giving rise to the optical flux (e.g. Idi Matteo et alj 1 19991 : 
ICorbel & Fended |2002| ; liioman et alj|2005l) , while the X-rays are 
thought to originate as a result of Compton ups cattering of seed 
photons either from the disc or from the jet JReis et alj 120081 : 
iMarkoff et al.ll2005[), or as a direc t result of optically-thin jet syn- 
chrotron (e.g. iMaitra et al.l 120090 . What new information do we 
have on the state of the source during our 2007 post-outburst 
low/hard state observations from the broad-band timing and spec- 
tral analysis above? 



Figure 17. Optical-only and X-ray-only CCFs. (Top) The plot shows the 
r' vs. g' CCF from Night 1 (light curve time resolution of 50 ms) as the 
black histogram, with lags refemng to the r' band. The orange histogram is 
the same CCF mirror-imaged about zero lag, i.e. lags referring to g' . Some 
asymmetry is apparent. (Bottom) This shows the 5-20 keV vs. 2-5 keV 
CCF from Night 1 - the same as in the top panel of Fig. ll5l but this time 
shown zoomed in and computed on a fine 2^** s resolution. The grey error 
bars denote the standard deviation amongst the CCFs computed from many 
segments. 



5.1 Spectral energy distribution model comparisons 

A blue optical spectrum and several emission lines are seen in op- 
tical spectroscopy (Fig. |3] §[3TT]l. The spectral slope is bluer than 
that expected from a canonical Ray leigh- Jean's tail , but is within 
the range observed in many other XRB systems JShahbaz et al.l 
[l996), where it is attributed to the outer disc. But the high source 
optical:X-ray flux ratio during the low/hard state that we probed 
does not support a disc-only scenario energetically. This is il- 
lustrated in Fig. [19] where the top panel shows t he absorption- 
corrected Swift spectrum of the source reported by iTomsick et al.l 
I2OO8, their spectrum 2 was quasi-simultaneous to our observa- 
tions). It is clear that the multi-colour disc (using the DISKBB im- 
plementation in XSPEC) vastly underestimates the flux of the opti- 
cal counterpart. Constraints on this disc component are weak, as the 
sensitivity of the Swift X-Ray Telescope (XRT) drops fast below 1 
keV, in addition to the unknown absorption. But these uncertainties 
cannot account for the large excess suggested for DISKBB (dashed 
line in the plot). The implied energy budget of a disc scaled up to 
the observed optical flux is about 50 per cent of the Eddington lu- 
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Figure 18. X-ray vs. optical multi-band cross-spectram components. The 
top plot shows the X-ray and optical (r'-band) PSDs and fits from Figs. 1101 
and 111! The following panels show the Coherence function, phase (</<), 
and time lags, respectively, going down. Red, green and blue symbols refer 
to the r', g' and u' bands respectively. Positive lags imply that optical is 
delayed. In the time lags plot, circles denote lags computed according to 
4>/2TTf. Empty circles with dashed eiTor bars denote negative lags, plotted 
at their absolute (positive) values on the log-log scale. For the three highest 
frequency bins, the triangles denote the time lags computed assuming that 
the phase lag increases smoothly above n (see text). Nights 1 and 2 have 
been averaged over the common frequency range for r' and g' . 



minosity for a 6 M© BH, completely inconsistent with the low/hard 
state. Instead, one may consider the DISKIR model JGierliiiski et al.l 
l2009l) . in which the disc is irradiated by central X-rays, and a frac- 
tion of the incident energy is thermalised to lower energies span- 
ning the ultraviolet to the infrared. Fig. [19] shows such a model 
(hereafter referred to as 'diskirl') with an inner disc temperature 
and high energy Comp t onisat ion tail parameters matching those 
found bv lXomsick et al.l ( 120081) . The main parameters are the frac- 
tion of reprocessing occurring in the outer disc (faut=0-l) and the 
Comptonised-to-disc luminosity fraction (Lc/Ld=lO). Other pa- 
rameters include the outer disc radius (rout=10'*'^rin) and a low 
fraction of reprocessing in the inner disc. But the problem with 
such a disc-only scenario is still the high optical:X-ray flux ratio, 
which requires a large fraction of reprocessing to occur in the outer 
disc. The above value of /out implies a very high efficiency of re- 
processing in the outer disc, orders of magnitude more t han in other 
systems (e.g. iGierliiiski et al.l2009l : lchiang et al.ll2010l ). Addition- 
ally, the spectral curvature of the outer disc drop-off curves in an 
opposite fashion to that of the broken power-law optical contin- 
uum (cf. Fig. [Sj. Finally, the very fast time variability, and anti- 
correlations (discussed further below) also point to the presence of 
other physics. 

Alternatively, pure jet models have recently been discussed in 
the litera ture. In Fig.|19l we show the results of one of the latest of 
these, bv lMaitra et alj ( |2009|), in which a bright hard state outburst 



SED from 2 002 reported by lHoman et alj ( l2005h is modelled (their 
'Model 2'). iMaitra et alj found that the source jet power, which 
could fully account for their observed optical flux, was higher than 
on any prior occasion when simultaneous broad-band observations 
were carried out. But the model comparison with our data in Fig. 1 191 
is seen to over-predict the observed X-ray flux by a factor of ~3 at 
10 keV, when normalised to the optical data (only a slight overall 
renormalisation by a factor of 0.4 has been applied to match to our 
2007 optical flux). In other words, the source optical:X-ray flux ra- 
tio in our observations is higher than that seen when the source was 
previously found to be fully jet-dominated. We note that no simul- 
taneous radio measurements of GX 339-4 around our observing 
dates have been p ublished (though a compact jet was detected; cf. 
iTomsick et alj|200 8), so we have used the radio/X-ray correlation 
of ICorbel et alj ( 120031) to extrapolate the 3-9 keV X-ray flux to 8.6 
GHz, and this monochromatic flux is also plotted in the figure and is 
found to agree with the normalised jet spectrum. Complete jet mod- 
elling is beyond the scope of the present paper, but the plot clearly 
shows that to be successful, it will have to incorporate a decreased 
contribution of the power-law tail electrons at higher energies. Ad- 
ditionally, a separate ionized gas component cool enough to pro- 
duce the observed optical recombination and fluorescence emission 
lines would be required. 

From the above comparisons, we find that 1) pure disc 
(re)emission is unlikely to explain the optical flux energetically, 
and 2) the optical:X-ray flux ratio is higher than found by previous 
jet modelling when the source was already optically very bright 
compared to its X-ray flux. 



5.2 The timing behaviour and its implications 

The detailed analysis of the optical vs. X-ray timing behaviour pre- 
sented herein is completely new for GX 339^. In particular, the 
Fourier lags indicate the presence of at least two components to the 
optical variability. 
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Figure 19. Absorption-corrected SED (cf. Fig. |4) compared witli vari- 
ou s models. {Ton) The solid blue curves show the disc and power-law 
of iTomsick et alj (2008), and the dotted curves are their blurred pexriv 
and gaussian line reflection components. The dashed blue curve is the disc 
scaled up by a factor of 250. The pink triple-dot-dashed curve is an ex- 
treme DISKIR model ('diskirl'), whereas the triangles denote a more typical 
DISKIR model (which we call 'diskir2'), taken from Gierlifiski et al. ( 2009) 
for the case of XTE JI817-330. (Bottom) Jet model ('2') from Maitra et alj 
120091) from radio to X-rays: the red dashed and the purple triple-dot-dashed 
curves represent post-shock (outer jet) and pre-shock (jet base) synchrotron 
respectively, and the dotted dark blue curve is the self-Comptonised pre- 
shock synchrotron. The red triple-dot-dashed curve is the sum of these. The 
plotted radio flux is a prediction based on the X-rays, assuming the low/hard 
state relation between the 3-9 keV integrated and the 8.6 GHz monochro- 
matic fluxes, found bviCorbel et alj I i2003i) . 



5.2.1 High Fourier frequencies 

Above ~0.2 Hz, phase lags (<^) increase almost linearly as a func- 
tion of Fourier frequency up to «2 Hz, then switch to a lower limit 
of ~— TT and continue increasing again. The corresponding time 
lags have a complex shape, with sharp rise from 0.4 to 0.7 Hz, fol- 
lowed by a gentle roll-over to about 3 Hz, and a quicker decrease 
thereafter. But the phase lag switch in this case is simply the result 
of the fact that it is defined within the range of -vr to tt (§ I4.4.2I >. 
For a smooth increase above the upper limit of tt, we can transform 
these negative (j> values at high Fourier frequency to an absolute 
phase value of 2-k+(J>. The recomputed time lags [i.e. (2n+(j))/2n f] 
then have similar values to the other high frequency time lags, and 
are plotted as the triangles in Fig. [TS] (bottom). The distribution 
suggests an approximately constant value of the time lags in this 



regime, close to the peak of the CCF time delay at 150 ms (Fig.llSll. 
Fitting a constant to the time lags above 0.4 Hz results in a com- 
patible lag of 0.145±0.004 s and 0.148±0.004 s for r' and g' re- 
spectively. The time lag uncertainties in this regime are very small, 
with the result that a simple constant is not statistically a good fit. 
But without more insight as to good physically-motivated models 
to try, we do not fit other detailed forms here. The important points 
are that: i) the range of variation of high Fourier frequency time lag 
values is only a factor of ~3 or less, which is small compared to 
the variation seen at lower Fourier frequencies; and ii) the time lag 
at the peak of the coherence function itself (at a Fourier frequency 
of 0.7 Hz) is measured to be 0.15±0.02 s and 0.14±0.02 s for r' 
and g' respectively, matching the peak of the CCF time delay. 

A constant time lag can be produced if a simple propagation 
delay operates between two signals. If one assumes that this cor- 
responds to a simple light-travel time, the distance between the X- 
ray and optical emitting regions is ~5000 Rg [= GM/c?} for 
A/=6 M0 (or a range of 2000-6000 Rg for the time lag variation 
range noted above). The narrow optical ACFs mean that the optical 
emitting region must be compact, ruling out the outer portions of 
the accretion disc as a possibility. Assuming that the optical fluctu- 
ations arise at the base of the jet (see Fig. |19t , the above distance 
is too large fo r typical jet accelera tion zones, which are thought to 
be <100 Rg JMarkoff et alJI 20051) . This association is much more 
reasonable if the fluctuations propagate instead at the sound speed 
(Cs) of the corona. In the case of a magnetized corona, this may be 
the Alfvenic velocity, which is expected to be only a fraction of the 
speed of light (c). For typical values of low/hard state plasma op- 
tical degths_Or~liasdrnagnetic compactness (Zs~0.5), Cs~0.01c 
(cf. iMalzac & Belmonl|[2009l , for constraints on Cyg X-1). In this 
case, the jet base lies at a distance of «50 Rg (with the time lag 
variation implying a range of 20-60 Rg)- On the other hand, the 
red variability rms (Table[3]and Fig. [8} may instead be arguing for 
the fluctuations originating in the optically-thin jet. In this case, the 
delay is a result of the distance between the optical and X-ray emis- 
sion regions traversed by fast relativistic jet plasma. The resultant 
lower limit on the jet elongation is then several thousand gravita- 
tional radii, as above. 



5.2.2 Low Fourier frequencies and QPO 

Below ~0.2 Hz, the lags again rise towards low Fourier frequen- 
cies, with the additional appearance of a 0.05 Hz QPO in the optical 
PSD associated with low coherence values. The long timescales at 
the lowest Fourier frequencies below O.OI Hz are probably associ- 
ated with weak reprocessing in the outer disc, as we discuss further 
in ij |5.4l Other contributions are also likely to exist over the interme- 
diate frequency interval. For instance, the phase lags of 7r/2 below 
0.1 Hz may contribute an anticorrelated or 'differen tial' response 
component (as suggested for XTE JI1I8-I-480; cf. iMalzac et alj 
I2003ll2004h . In addition, time lags of a few to ~20 s are too short 
for viscous effects to dominate, but may fit in with thermal insta- 
bilities at radii of several hundred Rg, although there is no obvious 
characteristic timescale to associate with the QPO. 

QPOs in the optical are relatively rare, so no general clas- 
sification exists. But how does this feature compare with the ex- 
tensive QPO detections in X-rays? The characteristic frequency of 
0.05 Hz, Q values of a few, and integrated rms values of ~2-3 
per cent (at least on the better Nights 1 and 2), suggest that this 
feature has properties similar to 'Type C low frequency QPOs 
classified bv lCasellae t al.' ('200 5J). The QPO frequen cy lies below 
the range studied bv lWiinands & van der Klij ( 119991) . where they 
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found a close relation between the low frequency broad-band break 
(vh) and QPO characteristic frequency. If our QPO agrees with this 
relation, then a break at fb ~0.005 Hz must be present. We note 
the presence of a weak upturn at the lowest frequencies in the op- 
tical PSDs in all filters (Fig. Illb consistent with this, though it is 
clearly important to sample lower Fourier frequencies with longer 
observations. 

A dynamical time ~20 s corresponding to the QPO charac- 
teristic frequency would mean a perturbation at a radius of ~2000 
Rg, which is too large for t he expected inner ra dius of the disc, 
even if it is recessed (cf. iTomsick et al.l 120081) . Epicyclic fre- 
quencies corresponding to Lense-Thirring precession dln gram et al.l 
[2009) could work. But one would then expect the QPO to appear 
in X-rays as well, unless it is somehow isolated from the coronal 
flow and strictly associated with the disc, to which RXTE is insen- 
sitive (although n ote that RXTE did d etect a low frequency QPO in 
XTE Jl 1 18-1-480: 'Hvnes et al.ir2003ii) . We do not detect an X-ray 
QPO, but have already noted in ii l4.2.3l that we cannot rule out the 
presence of an X-ray QPO identical to the optical one. On the other 
hand, the optical QPO carries ~20-25 per cent of the total source 
rms variability detected (cf rop?Lai=0-15-0.12; Tabled, whereas 
our derived limit of ^'x-rav^^-^^^ "-"^ ^"^y X-ray QPO implies an 
upper-limit of only 16 per cent of the corresponding X-ray source 
rms. Thus, the process generating any QPO common to both ener- 
gies must contribute a larger fraction to the optical variable flux, as 
compared to X-rays. 

With regard to other simultaneous optical and X-ray observa- 
tions of the source, no X-ray QPO detection has bee n reported at 
X-ray flux levels similar to ours. The timing study of lNowak et al.l 
( 11999 1) carried out at a flux level only two times higher than ours 
(§ l4.2.4t should have been much more sensitive because of the fact 
that all five PCU detectors on the PCA were operational. Although 
no targeted search at characteristic frequencies around ~0.05 Hz 
was performed, there is no obvious QPO in their presented fig- 
ures. On the other hand, a prominent low frequency X-ray QPO 
was found by Motch et al. (1983) using the Ariel 6 satellite in con- 
junction with their optical QPO, at a low/hard state X-ray flux level 
~5 times higher than ours. The X-ray QPO amplitude was reported 
to compare well with that seen in the optical. So clearly, the X-ray 
QPO properties are changeable even within the low/hard state. The 
LAXPC instrument on board the upcoming Astrosat satellite has an 
effective area of about 3 times that of the RXTEfPCA at hard X-ray 
energies ( lAgrawalll2006l) and will be an excellent tool for identify- 
ing any weak low-frequency QPOs during X-ray faint states. 



5.2.3 Overall coherence and PSD change patterns 

The high Fourier frequencies show a significant optical vs. X-ray 
coherence, mainly associated with the !-'max~0.9 Hz Lorentzian 
component (top panel of Fig. 1 181 and Table [Sjl- A slight bump be- 
tween 2 and 3 Hz also seems to coincide with the highest frequency 
Lorentzian above 3 Hz. At the very highest frequencies, the lack of 
significant X-ray power (above PSD white noise) is responsible in 
part for the loss of coherence. At lower Fourier frequencies, the 
cross-spectrum behaviour is complicated. The coherence decreases 
sharply from «0.3 to 0.1 Hz and the phase lags in this regime also 
undergo a sign reversal, in the sense that the optical leads the X- 
rays. The uncertainties are smaller than the absolute phase values. 
These effects are seen to coincide with a change in the dominant 
PSD variability components, with the coherence decline associated 
with the decrease of the 0.9 Hz optical Lorentzian. The negative 



lags occur in conjunction with the rise of the strongest X-ray PSD 
component around 0.1 Hz, and also the emergence of the lowest 
frequency zero-centred optical component (i'inax~0.06 Hz). Such 
a changeo ver in fluctuations can naturally cause a loss of coherence 
(e.g. Vau ghan & Nowaklll997h . especially if the origin of the low 
and high freque ncy Lorentzians is indepen dent, which is likely to 
be the case (e.g. I Wilkinson & Uttlevll 20091) . 



5.3 RMS spectrum and broad-band constraints 

What does the energy spectrum of variations look like? This can 
be examined by combining the results from the spectral and tim- 
ing analysis as follows. One first computes the fractional vari- 
ance by integrating the noise-subtracted X-ray and optical PSDs 
over any desired Fourier frequency range (from /i to /2). The ob- 
served fluxes (-Fband) in all bands are then multiplied by the respec- 
tive fractional rms values (i.e. the square-roots of the variances) 
to give the broad b and rms spectral energy distribution (see, e.g., 
iHvnes et al]|2003ij) : 



-f^rms — -t^h. 




P{f)df 



(8) 



where P represents the noise-subtracted PSD model fits for each 
band in rms-squared normalised units. We choose the Fourier fre- 
quency range of 0.4-4 Hz for this exercise because this is associ- 
ated with a relatively uniform optical/X-ray coherence (Fig. I18t . 
This range is dominated by the strongest optical PSD Lorentzian 
component, and in X-rays, includes the peak of the high frequency 
Lorentzian. Fig.l20lshows the absolute rms fluxes for both the fast 
r' and g' light curves, and for the 2-5 keV and 5-20 keV ranges in 
X-rays. The respective mean fractional rms values were measured 
to be 0.11, 0.09, 0.23±0.02 and 0.24±0.03, where the quoted er- 
rors are statistical la uncertainties determined from Monte Carlo 
sampling of the PSD fits to obtain an ensemble of rms measure- 
ments. These errors are very small in the optical, where systematic 
dereddening uncertainties instead completely dominate (overplot- 
ted in the figure). The dotted line is a power-law fit to the four rms 
points, and yields a power-law index of 0. 17±0.04 in the plot units 
of the figure (x^/dof=3/2), or an index of -0.83 in the commonly- 
used spectral energy density (= F^) units. Instead of the 0.4^ 
Hz Fourier frequency range, if we use the full PSD range from 
Figs.[TO]and[TT]for this exercise, we obtain a broad-band slope of - 
0.78±0.04. These values are significantly steeper than the slope in- 
ferred from spectral fitting to the full X-ray spectrum (ij |3.2t . where 
we found an energy index of l-r=-0.63(±0.01) or -0.66 (±0.04), 
depending on the model used. They are also much steeper than the 
slopes found in the optical (energy in dices of to - l-l; JJI3.U . 

A power-law was also found by Hvnes et al.l ( l2003al see also 
iHvnes et al.l200Q) for the rms spectrum of XTE J 1 1 18-1-480, and as- 
cribed to a single optically-thin jet synchrotron component extend- 
ing from the near-IR to X-rays. Our findings suggest that if there 
is a such a single broad-band varying optically-thin component in 
GX 339-^, it must have a spectrum different to that describing the 
mean emission in both optical and in X-rays. Additional compo- 
nents with different spectral and variability characteristics, such a 
jet base, magnetic coronae or inner hot accretion flows in the optical 
(e.g. Merloni et al. 2000 ; Yuan et aL 2005 ; Maitra et a l. 2009) , and 
Compton upscattering of jet and disc seed photons in X-rays (e.g. 
iMarkoff et al.ll2005l ; lMakishima et al.ll2008h . may then account for 
the remaining spectrum. 
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But the problem with such an interpretation is the low abso- 
lute value of the optical:X-ray coherence of «0. 1, which means 
that a single component is not necessarily a good representation of 
the broad-band variability. Were such a component to dominate the 
variability, an optical:X-ray coherence of close to unity would be 
expected. Instead, it is to be noted that the rms values in both X-ray 
bands are almost identical, meaning that the X-ray-only rms spec- 
tral slope is close to the mean flux X-ray spectrum, and is flatter 
than the single broad-band rms spectral slope inferred above. In the 
optical, the g' rms is lower than that in r', indicative of the red rms 
slope discussed in § 14.11 and Fig. [8j still the r'~g' slope is again 
significantly flatter than the broad-band one. The dotted lines in 
Fig.|20]denote the X-ray-only and optical-only rms energy slopes. 
These flat slopes support the scenario whereby separate physical 
components dominate the variability in the optical and in the X-ray 
bands, with the positive CCF implying some underlying cormection 
between the two. Note that this inference holds true in the optical 
despite the large dereddening uncertainties, i.e. for all Ay values 
ranging over 3^ mags at least. 

In X-rays, the variable component describes the mean flux 
spectrum well. This may be powered either by an optically-thin 
jet, or Compton upscattering. We cannot distinguish between these, 
but we note that a disc extending down to an inner radius to ~10 
Rg is seen to be pr esent in soft X-rays below 1 keV by Swift 
JTomsick et al.l 12008 1. This, as well as the prominent jet base or 
inner coronal flow which likely dominates the optical (see below), 
can all easily serve as seed photon sources for Comptonisation (cf. 
iMakishima et alj2008l for a similar discussion with regard to Cyg 
X-1). 

In the optical, a straightforward cause of the lower g' rms vari- 
ability as compared to r' may be a steady component such as an 
irradiated disc, which is expected to emerge mainly at blue wave- 
lengths. An estimate of the disc's maximal contribution at g' can 
then be simply computed by using the excess r' rms (over g'), and 
turns out to be «20 per cent. This is likely to be an upper limit, 
because other possible components such as a jet base, magnetic 
coronae and a i r mer hot acc retion flow (e.g. IMerloni et al.ll200d : 
lEsin et al.ll200ll ; lYuan et alj|2005,) which may be present also peak 
in the blue. DISKIR parameters more typical than those required by 
the extreme diskirl model discussed in ij 15.1 l ean easily produce a 
component consistent with this. In Fig . [19| we overplot a m odel 
from the high end of the range found bv lGierlinski et al.l ( 120091) for 
the source XTE J18 17-330, with the fraction of reprocessing on the 
outer disc /out=0.01 and the irradiation fraction at the inner disc of 
/in=0. 1 (model hereafter referred to as 'diskir2'). This results in an 
optical disc:total flux ratio ~10 per cent which is also consistent 
with a canonical disk heating scenario discussed i n next section, 
and sa tisfies the Swift soft X-ray disc constraints of iTomsick et al.l 
( 120081) . Additionally, in this scenario, the broken optical spectral 
continuum detected by FORS may be explained by the appearance 
of the disc at high frequencies, for which the power law fits of Fig.|3] 
suggest a contribution of ~16 per cent at the bluest wavelengths of 
4000 A. 

Further testing for the existence of any additional underly- 
ing single broad-band component will requir e extending such an 
rms analysis to low energies in the near-IR (cf. lHvnes et al.ll2003al : 
ICasella et alj|201(]|) , because this regime is not expected to have 
strong contributions from the disc or inner hot flow and jet base 
discussed above. 

That the optical and X-ray CCFs show a positive peak means that 
the variability in the bands is connected. At the same time, the rms 
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Figure 20. Optical and X-ray absolute 0.4^ Hz rms spectral fluxes (lower 
four points) compared to the total observed spectra. The fractional rms val- 
ues derived from the PSDs in the r' (centred at the effective wavelength of 
6250 A), g' (4800 A), 2-5 keV (3.5 keV) and 5-20 keV (10 keV) bands, 
from left to right, are 0.11, 0.09, 0.23±0.02 and 0.24±0.03, respectively. 
Systematic dereddening uncertainties to the optical are shown in green. The 
pink dashed line is a fit to these four points and has a slope of 0. 17±0.04 in 
the plotted log-log units. The orange and blue dotted hnes connect the two 
optical and two X-ray rms points respectively; the slopes of these two are 
flatter than that of the pink line. 



analysis in this section suggests that distinct physical components 
dominate the variability in each band. A self-consistent scenario 
that satisfies the main constraints of our work may then be as fol- 
lows. 

X-ray and optical studies have now shown that perturbations 
in the outer parts of the accretion flow propagate towards the centre, 
seeding and coupling fluctuations in all inner re gions, and resultin; 
in a linear relation between rms and flux (e.g. luttlev et alJuOO: 



lGandhill200ft , and references therein). If the disc (required by the 
Swift data) is magnetically-threaded and anchors a jet (which could 
explain the bulk of the optical flux), then perturbations will travel 
along the field lines leading to line reconnection accompanied by 
X-ray flares in the innermost regions. These will be Comptonised to 
high X-ray energies with negligible time delay to produce the hard 
X-ray power-law and high X-ray-only coherence values that we 
observe. Upon reaching the jet, the perturbations are expected to 
trigger (cyclo-)synchrotron emission, leading to the optical emis- 
sion and hence the correlated optical :X-ray CCFs that we detect. 
At the same time, the slower and broader CCF anti-correlation 
humps can result from decreased coronal cyclo-synchrotron caused 
by magnetic field dissipation on the X-ray flaring timescales (see 
further discussion in Paper I). Thus, such a propagation sequence 
can generate the CCF features that we observe. 

In this picture, the low peak values of the coherence function 
and CCFs are a result of interaction between at least two distinct 
components: the jet (generating the optical power) and the corona 
(emitting X-rays). It may be that propagation and energy transfer 
events between the two components are episodic, or that the trans- 
fer occurs via some 'filter' (e.g. magnetic reconnection, or change 
in field line dynamics) which ends up diminishing coherence (see 
also Paper I; Durant et al. 2010 in prep.). Propagation of pertur- 
bations will also naturally lead to linear rms-flux relations in both 
X-rays and in optical. The perturbations may arise directly within 
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the hot flow itself (e.g. iDone et al.l2007h as long as the flow is large 
enough to encompass the broad range of timescales required by the 
rms-flux relation. 

If some form of cyclosynchrotron is indeed important for the 
optical variability (whether in a pre-shock jet or a corona), one can 
estimate a minimal magnetic field (B) energy density from one of 
the observed prominent flaring episodes, assuming that the B field 
is in equip artition with the radiati on field. For this estimation, we 
use Eq. 1 of ldi Matteo et alj ( ll997h with the following assumptions: 
1) dissipation of energy within the disc is negligible compared to 
that in the jet base or corona, as should be true for the low/hard 
state; 2) a single emission region is responsible for powering the 
flare, and can be described as a uniform electron cloud with di- 
ameter ~100 milli-light-seconds, equal to the fast decay times that 
we observed (Fig.|7J; and 3) an integrated radiative peak flare out- 
put of ~10^'^ erg s"^ RiO.OOlxLEdd- This yields a B field den- 
sity of 5 X 10'^ G, which is likely to be a very conservative lower 
limit, given that much faster optical flaring timescales have been 
detected previously qMotch et al.ll 19820 . and that the corona may 
be patchy. Both these effects will considerably decrease the emit- 
ting region size, and in turn increase the B field density. Hence, the 
observed optical flaring radiative losses are easily consistent with 
strong field energy dissipation events (e.g. War dzinsk i & ZdziarskI 



1 



l200d : ldi Matteo et alJl997l : lFabian et anil982l) 



5.4 On the reprocessed component 

J. 4. 1 The fast variability 

The PSD plots of Fig.[9]show that there is consistently higher power 
in the X-ray variability as compared to the optical. So could the 
optical variability be a result of reprocessing of faster X-ray vari- 
ations at frequencies above the optical PSD peak (or even higher 
than those probed by us due to white noise limitations)? In Paper I, 
we argued against reprocessing dominating the rapid optical vari- 
ability on times of <: 1 s, at least as described by simple linear 
transfer functions. The reasons are several. 

Most obviously, the anti-correlation between optical and X- 
rays is the opposite of the expectation under a reprocessing sce- 
nario, unless on e is willing to invoke negative transfer functions 
JHvnes et alT2003a) . 

Then there is the very short time lag and narrowness of the 
peak of the positive CCF component for the fast light curves. For 
an orbital period of 1.7 days, the corresponding binary separation 
is ~25 light-seconds. Our observations on the three alternate nights 
were carried out with an inter-observation spacing of about 47 and 
49 hrs, respectively (see Fig.[2l(, or relative phases of 0.15 and 0.38. 
For X-rays irradiating the companion star followed by an optical 
light echo, the expected time lags are then much longer than the 150 
ms peak delay observed on all the nights. In any case, GX 339-4 is 
thought to be a low inclination binary which increases the minimum 
expec ted time delay to a(l-sim)=19 s for i=15 degrees jWu et al.l 
1 200 I I). Reprocessing in the outer disc, instead of the companion 
star, ought to peak on much longer times and also be smeared out, 
and so is clearly unrelated to this fast variability. 

Another clue is the strength of the X-ray vs. optical CCF on 
the shortest lags, which is strongest (weakest) in the r' (u') filter, 
the opposite of the expectation based on disc reprocessing. 

Finally, the optical ACFs are wider than the X-ray ones, in 
all optical filters (see discussion in § 14.31 ). This is another way of 
saying that the peak of the optical variability power is skewed to- 
wards higher frequencies as compared to the X-ray peak. The PSD 



decomposition in § |4] showed a dominant optical component with 
i^max~l Hz, whereas the X-ray PSD peaks around 0.1 Hz. 

None of these arguments, by themselves, rule out a more com- 
plex reprocessing origin for the fast optical peak delay, especially 
given that the low X-ray source count rates limit our ability to de- 
tect high frequency variability power. But such an origin also re- 
quires that the reprocessor be 1) compact; 2) stable over several 
days at least; and, 3) able to isolate the fast X-ray variable compo- 
nent from the slow one. This can be discerned directly by eye from 
Figs. |5] and (6) where the X-ray light curves have very significant 
broad flaring components during which the average source flux is 
increased by factors of a few, over times of several seconds. The 
simultaneous optical section shows no obvious positive response to 
this. If anything, the optical shows local minima over a period of 
several seconds around the X-ray flares in both figures. So if part 
of the X-rays are being reprocessed on sub-second timescales on 
some reprocessing region, the slower and broader flares are not si- 
multaneously seeing the same reprocessing region. This scenario is 
inconsistent with canonical disc reprocessing models, at least. 



5.4.2 The slow variability 

Disc reprocessing is expected to be stronger at bluer wavelengths, 
because the hot inner regions (of the outer disc) can intercept a 
larger fraction of incident photons from a central X-ray source. So 
bluer filters ought to show stronger signal on timescales commen- 
surate with the outer disc. This is consistent with the behaviour 
of our slow CCFs in Fig. [TS] and also the high u' coherence val- 
ues and long time lags of Fig. [TS] So it is natural to associate 
this slow component with disc reprocessing. Such a component 
increasing towards the ultraviolet has also been ide ntified in the 
case XTE J 1 1 1 8+480 bv lHvnes et Zl J2003al . see also lMalzac etall 
[20041). The Bowen blend that we observe in the optical spectrum is 
also thought to be a result of reprocessing (fluorescence), though 
its location, whether in a disc, or wind, or a low density gaseous 
component, is unclear (cf. ISchachteret alJI 19891 ; IWu et al.ll200ll) . 
Echo tomography of the line during subseq uent active low/hard 
states should help to set tle the issue (e.g. lO'Brien et al.l 120021 ; 
iMuiioz-Darias et alj2007t) . It is worth mentioning that if the Bowen 
blend arises in the disc, our limits on the optical disc contribution 
(i) l5.3t imply that the equivalent width of the blend against the disc 
continuum must be at least 5 times larger than the observed equiva- 
lent width. Finally, we note that exactly which Fourier frequencies 
correspond to reprocessing is also unclear because of the changing 
phase lags below 0.1 Hz (see Fig. |18t; only the very low Fourier 
frequency bin shows an optical phase lag of close to zero with re- 
spect to X-rays. 

Energetically, the reprocessed (disc) component must be weak 
unless the reprocessor is highly atypical, as already discussed in 
i^ 15. II with regard to the SED models of Fig. [19] Under canonical 
disc heating scenarios, what is th e average optical flux expected? 
Ivan Paradiis & McClintockl (19940 have shown that in'adiation by a 
central X-ray source on the outer portions of an optically thick disc, 
where visual surface brightness varies as T^ (T being the local disc 
temperature), leads to the optical luminosity varying in proportion 
to E=(Lx/I'Edd)^ ^^h , where P is the orbital period in hours. 
For our observations of GX 339-4 carried out at a low Edding- 
ton rate, E~l. For this value of S, an absolute visual magnitude 
of Mv=l.(> is predicted as a result of reprocessing. Our inferred 
Mv=-1.02 (§[33, an excess flux of a factor of 11. The implied 
reprocessed flux is close to that of the irradiated disc diskir2 model 
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shown in Fig. [191 Dereddening uncertainties only change this by 
a factor of two, while distances of 5-15 kpc imply a range for the 
excess flux factors (above the reprocessing prediction) of 7-21. 



5.4.3 Lags longer than ^ 10 s? 

It is clear from the discussion so far that at least two timescales, 
associated with the fast (~150 ms) and the slow (~10 s) variabil- 
ity lags respectively, and one more related to the QPO character- 
istic time (~20 s), are important in GX 339^. What about other, 
even longer timescales? One way to highlight these may be to filter 
the light curves on timescales better matched to the t ime delay be- 
ing so ught. We did this by following the procedure of lMalzac et al.l 
( 120031) . In short, the CCFs were computed after keeping only the 
time-scales in a specific range of frequencies. High frequency noise 
was removed by smoothing with a box car filter. Low frequency 
noise was removed by dividing the light curve by a piecewise lin- 
ear trend, i.e. a linear interpolation of the light curve on the longest 
time-scale retained. The results are displayed in Fig. 12 1 1 

This exercise highlights various features seen in the raw (un- 
filtered) CCF shown at the bottom. In particular, the main positive 
peak appears at longer times, matching the expected delay on sev- 
eral seconds expected from reprocessing. Some wiggles due to the 
20 s QPO are also apparent. On the other hand, the negative correla- 
tion components are also seen to shift smoothly to longer (absolute) 
time delays, whic h is not expected. Su ch 'self-similar' behaviour 
was also noted bv iMalzac et^ ( l2003h for XTE Jl 118+480, and 
used a basis for the stochastic fluctuation model of magnetic 
reservoir energy release which successf ully explains the sp ectro- 
temporal characteristics of that source dMalzac et al J 1 20040 . But, 
unlike XTE Jl 118+480 where afl time lags from 0.01-10 s con- 
tri buted smoothly ove r the full Fourier frequency range (e.g. Fig. 8 
of 'Malza c et al.ll2003 '). the main contributions in GX 339-4 are the 
150 ms, and the ~10 s component (bottom panel of Fig. 118b. No 
other important features are apparent on timescales of ~ 100 s, at 
least. 



5.5 Comparison with other sources 

Several distinct components (including both reprocessed and di- 
rect variabilit y) contribute to t he optical power of XRBs in 
general (e.g. iRussell et al.l 120061) . There are indications of non- 
reprocessed optical and IR (OIR) variability i n several sources now 
teikenberry et al.' 1998; Uemura et al. 2003; iKanbach et ai]|200ll ; 
[Purant et al...2008. ; .Gandhi et al...2008.) . GX 339-4 has one of the 
highest fractional varia bility amplitud es and strongest non-linear 
flaring amongst these JGandhil 120091) . Fast variability was typi- 
cally found in the low/hard state in all sources, though OIR ob- 
servational coverage in other states still remains very sparse on 
rapid timescales. In all cases, the sources exhibited bright coun- 
terparts, and rather high OIR:X-ray flux ratios. Which low/hard 
state property drives this non-reprocessed component is currently 
unclear. The jet seems to play a key role in both XTE Jl 1 18+480 
and GX 339-4, but SWIFT J1753.5-0127 possesses only a weak 
jet, with the corona and inner disc playing more important roles 
JDurant et ai]|2009l) . On the other hand, the optical ACFs of both 
the former sources were narrower than the corresponding X-ray 
ones, which is not true for SWIFT J1753. 5-0127. This does suggest 
that the very fastest flaring may be associated with the jet. Finally, 
the possibility of jets with strong magnetic fields and correspond- 
ing OIR synchrotron emission, but with a quenched radio emission 
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Figure 21. X-ray vs. r' CCFs filtered in order to highlight various 
timescales. Only times within the range labelled at the top left were retained 
in each case. Times shorter than the smallest labelled ones were smoothed 
out using box car filtering. Those longer than the longest labelled one were 
filtered by piecewise Unear normalisation. Note that smoothing can artifi- 
cially increase the absolute correlation strength on the y-axis. 
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(hence makin g their detection less se cure observationally), cannot 
be discounted l iCasella & Pe'en2009l) . It may also be the case that a 
direct OIR variable component is generally present in other XRBs 
as well, but is smothered by the typically brighter reprocessed com- 
ponent. Observations in the high/soft state sensitive enough to be 
probe beneath the expected reprocessing level are needed to con- 
firm this. 

XTE Jl 118-1-480 is the only other source with detailed rapid 
optical vs. X-ray coherence and lag analyses published so far As 
mentioned in the previous section, one key distinction between 
GX 339-4 and XTE 111 18-1-480 is that in the latter case, a con- 
tinuous range of time lags contribute to the coherence function, 
whereas we find only a fast (~150 ms), and a slow (~10 s) 
timescale to be important for GX 339-4. In terms of their physical 
characteristics, the most important difference between the sources 
is their orbital periods: GX 339-4 is a long period binary with 
P =1.75 d whereas XTE Jl 1 18-1-480 has a short period of only 0.17 
d JMcClin tock et al. 2001), which implies a correspondingly larger 
disc extent for the former, assuming that the Roche lobe is filled. 
Whether or not a larger disc can somehow suppress the contribu- 
tion from a continuous range of timescales in the case of GX 339-4 
remains to be seen. One way to test this would be to compute the 
coherence and lags for other systems: e.g. SWIFT J1 753. 5-0127 
has th e shortest period currently-known for any XRB JZurita et al.l 
l2008l) . 

Finally, the detection of a low frequency QPO only in the opti- 
cal and not in X-rays (to within the limits discussed in ij |5.2.2t may 
be unique to GX 339-4. XTE Jl 118-1-480 sh owed a QPO in both 
optical and X-rays during its 2000 outburst dHvnes et alj|2003al 
Spruit et al. 2010 in prep.). That the QPO in GX 339^ may be 
more stabl e in the optical, as o pposed to X-rays, has also been dis- 
cussed by iMotch et alj ([1985). This argues against an association 
with some instability at a radius characterising the edge of an inner 
hot flow, as these would produce X-ray QPOs as well. Instead, the 
QPO origin may be related exclusively to the optical emission re- 
gion. If this is the jet, then the characteristic time (and hence size) 
scale of the QPO may be telling us about the extent of the region 
where perturbations are being pumped into the jet. In this case, 
searches for similar oscillations at radio frequencies with future 
sensitive telescopes (e.g. the E VLA) could be a good discrimina- 
tor We note that Po olev & Fen der (1997) found radio QPOs in the 
strong jetted source GRS 1915-1-105 on timescales associated with 
IR flares in this system. On the other hand, GRS 1915-(-105 also ex- 
hibited soft X-ray oscillations on similar times, so the mechanism 
is not unique to the radio in this case. 



6 SUMMARY 

Our main results may be summarised as follows: 

(i) We have carried out rapid timing observations of GX 339-4 
simultaneously in optical and X-rays, on a best optical time resolu- 
tion of 50 ms. These are the fastest such observations of the source 
in the optically-dim (V~17) low/hard X-ray state. 

(ii) Simultaneous optical spectroscopy shows the presence of a 
blue, broken power-law continuum, and a host of recombination 
lines as well as the Bowen fluorescence blend (§[33J- 

(iii) The light curves reveal the presence of fluctuations on a 
range of timescales in both optical and X-rays (Fig.[5](. The fastest 
optical flares that we detect have decay times as short as ~100 ms, 
while the strongest flares can increase the observed flux to more 



than twice the mean level on these short times (Fig.|7]l. The multi- 
filter rms values suggest that flaring is redder than the absolute flux 
spectrum (Fig. [8). 

(iv) Simultaneous multi-filter power spectrum fits extending 
over three decades in Fourier frequency, and to well beyond 1 Hz 
in g' and r' (and slower in u') are presented for the first time, and 
show the presence of a low-frequency QPO in all filters (§ 14.21 
Figs. [ToUOt . No corresponding sharp feature carrying a similar 
fraction of the total X-ray variability is seen in the simultaneous 
X-ray data (ii l5.2.2t . Most of the optical variability power emerges 
around ~ 1 Hz. The X-ray rms is larger than in the optical, but peaks 
at lower frequencies around ~0. 1 Hz (Fig.|9]l. 

(v) Auto and cross-correlation analyses shows consistent results 
over several nights and in all filters, with the optical ACFs being 
narrower than the white-noise-corrected X-ray ones (Fig.ll4t. The 
CCFs in r' and in g' show a sharp and fast delayed peak with re- 
spect to X-rays (Fig. I15t . The u' data is slower, but also shows 
evidence of the presence of fast fluctuations (Fig.ll6t. 

(vi) Other than XTE Jl 1 18-1-480, our work presents the first de- 
tailed analysis of the fast optical vs. X-ray time lags and coher- 
ence function in an XRB (§ 14.4. 2t . The clearest feature of the 
coherence function between optical and X-rays is a broad peak 
around a Fourier frequency of 1 Hz, and this is associated with a 
time lag that is approximately constant at ~150 ms, the peak CCF 
delay. Changes in coherence values appear to be associated with 
changeovers in the dominant optical (or X-ray) PSD Lorentzian 
components with Fourier frequency (Fig.ll8t. 

(vii) We find a large optical-to-X-ray flux ratio (§[33J which is 
apparently higher than on any previous occasion when simultane- 
ous broad-band data was observed. This cannot be easily explained 
by a disc unless it is reprocessing an atypically-high fraction of in- 
cident X-rays. Furthermore, it requires that any jet model of the 
source in this state produce a high fraction of pre-shock optical 
synchrotron in the base as compared to the contribution from non- 
thermal optically-thin X-rays (§[5TT]l. 

(viii) Several lines of reasoning are presented against reprocess- 
ing of X-rays dominating the optical power and the rapid variabil- 
ity, at least according to canonical linear transfer function scenarios 
in the outer disc. These include: the low source X-ray Eddington 
fraction; the sharpness of the fast CCF peak and its small time lag; 
the anti-correlation seen in the CCF on times of <few seconds; 
and, the narrowness of the optical ACFs with respect to X-rays. 
Any other complex reprocessing scenario, if true, must incorporate 
atypically strong reprocessing ratios, and must also be able to iso- 
late the fast X-ray flares from the slower ones (§|5ATJ. 

(ix) Standard heating of the outer disc predicts a weak repro- 
cessed flux about 10 times lower than the total observed power. 
The slower (~10 s) variability component highlighted in the CCFs 
computed from the heavily binned light curves, and also in the large 
coherence values and lags at low Fourier frequencies, increases 
at blue (u') wavelengths, and is consistent with disc reprocessing 
(§|5A2l(. 

(x) The rms energy spectrum of fast fluctuations is best de- 
scribed by an X-ray variable component matching the overall en- 
ergy spectrum, and a separate optical component redder than the 
total optical energy spectrum (§ 15.3b . The contribution of a steady 
component such as a disc is constrained to be <20 per cent in g', 
and a typical DISKIR model extending from soft X-rays can easily 
satisfy this. A model in which perturbations related to accretion in- 
stabilities in the outer disc propagate to the inner flow and then to 
the jet, triggering hard X-ray flares via Comptonisation and optical 
cyclosynchrotron emission, can qualitatively explain the main re- 
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suits of our analysis, and can also produce linear rms~flux relations 
in both bands. The low optical :X-ray coherence peak value of ~0. 1 
is a result of interaction between these independent components. 
The B field necessary to produce the fast and strong optical flares 
that we observe must be >10* G. Furthermore, in such a scenario, 
the 150 ms time lag is related to the propagation delay between the 
X-ray corona and the jet optical emission regions. 



7 CONCLUSIONS 

Our wide optical/X-ray spectral and timing simultaneous coverage 
gives new insight on the processes occurring on fast timescales in 
the inner accretion regions of this important source. Fast optical 
timing is beginning to probe interactions of the jet, corona and disc. 
This provides good starting material for detailed modelling of the 
accretion environment. But many questions remain to be answered, 
including the origin of the optical QPO, the correlation of the ob- 
served lags with the various variability components found in the 
PSDs, the origin of the optical emission lines, and the underlying 
driver of the fast fluctuations and observed power in GX 339-4, 
and also other sources. The present inferences are based on rel- 
atively short datasets, with the best simultaneous dataset having 
a total duration of less than one hour. Longer follow-up observa- 
tions in various states of the sources should enhance the constraints 
presented herein. Indeed, several systems are now seen to show 
interesting variability characteristics on fast times, and such work 
ought to continue with new instruments being mounted on large 
telescopes and offering such modes on an increasing basis: e.g. 
FORS/VLT high time resolution mode, S ALT/S ALTICAM and Ul- 
traspec, among others. Further insight could come from improving 
the time resolution of the observations, or through detailed dynam- 
ical cross-correlation analyses utilising each segment of the light 
curves separately (cf. Durant et al. 2010, in prep.). Coordination 
with the infrared offers another excellent opportunity, as this may 
be abl e to probe optically-thick jet emission. Indeed, Casella et al.l 
( 120101) report the detection of fast infrared flaring in GX 339-4 
with the infrared also delayed by ~100 ms with respect to X-rays. 
Finally, instruments such as the LAXPC on board Astiosat will pro- 
vide increased hard X-ray sensitivity as compared to RXTEfPCA 
JAgr awal 2006) and thus also help to break some of the model de- 
generacies discussed herein. 
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Band 



Full PCA 



2-5 keV 



5-20 keV 



Lorentzian 
component 



(Hz) 



A 
(Hz) 



(Hz) 



X^/dof 



0.16(3) 
0.49(1) 
0.34(3) 

0.18(10) 
0.41(11) 
0.37(2) 



0.012(2) 



0.6(5) 

0.010(5) 

0.015(51) 





0.14(3) 0.012(2) 
0.51(1) 

0.38(6) 0.2(8) 



0.011(4) 

0.113(9) 

2.1(2) 

0.015(12) 
0.10(1) 
1.7(4) 

0.009(4) 
0.114(9) 

2.4(3) 



0.017(3) 

0.113(9) 

2.1(2) 

0.018(10) 
0.10(2) 
1.7(4) 

0.015(3) 

0.114(9) 

2.4(3) 



0.6(2) 



0.14(11) 

0.3(3) 

0.1(2) 





37/48 



53/48 



0.7(3) 
46/48 

0.04(17) 



Table 4. Multiple Lorentzian fits to the X-ray PCA PSD averaged for all three nights, listed in order of increasing i-'max 
P{iy) = r^A/7r/fA^ + (i/-i/n)^1. as defined bv lBellonieralH2002t) .i/^, 



The Lorentzian functional form is 
uo/2/A, where r, i^max and Q refer to the integrated 



= ^!/2 + A2 and Q 

rms (over the full range of -oo to +00), characteristic frequency, and quality factor respectively, n numbers in brackets denote the Icr uncertainties on the last 
n digits. 



Band 


Lorentzian 


r 


UQ 


A 


l^inax 


Q 


X^/dof 




component 




(Hz) 


(Hz) 


(Hz) 






r' 


1 

2 
3 
4 


0.029(4) 
0.074(7) 
0.192(3) 
0.055(8) 


Night 

0.049(1) 





2.9(3) 


i 

0.004(2) 
0.06(1) 
0.90(6) 
1.8(2) 


0.049(1) 
0.06(1) 
0.91(6) 

3.4(3) 


6.2(2.9) 





0.8(1) 


95/82 


9' 


1 

2 
3 
4 


0.024(4) 
0.058(6) 
0.154(3) 
0.057(6) 


0.049(1) 





2.9(3) 


0.004(2) 
0.06(1) 
0.89(6) 

2.2(2) 


0.049(1) 
0.06(1) 
0.89(6) 
3.6(3) 


5.9(2.6) 





0.7(1) 


94/82 


u' 


1 
2 
3 


0.05(1) 
0.020(4) 



0.050(2) 


0.02(1) 
0.005(2) 


0.02(1) 

0.050(3) 

>0.3 




5.5(2.7) 


22/19 








Night 2 








r' 


1 
2 

3 
4 


0.055(6) 
0.032(3) 
0.132(10) 
0.177(7) 




0.053(1) 






0.011(6) 
0.004(1) 
0.31(5) 
1.31(10) 


0.011(6) 

0.053(1) 

0.31(5) 

1.3(1) 




6.1(1.8) 






76/67 


9' 


1 
2 
3 
4 


0.052(6) 
0.026(3) 
0.076(22) 
0.155(6) 



0.053(1) 
0.08(9) 




0.016(7) 
0.005(2) 
0.22(7) 
1.23(9) 


0.016(7) 
0.053(1) 
0.23(7) 
1.23(9) 




5.1(1.6) 

0.2(2) 




62/67 








Night 3 








r' 


1 
2 

3 


0.054(4) 
0.178(21) 
0.141(26) 


0.052(2) 




0.011(3) 
0.47(8) 
1.5(4) 


0.053(2) 
0.47(8) 
1.5(4) 


2.3(5) 




43/38 


9' 


1 

2 
3 


0.036(5) 
0.077(9) 
0.160(4) 


0.052(1) 




0.006(3) 
0.13(4) 
0.86(6) 


0.053(2) 
0.13(4) 
0.86(6) 


4.2(2.1) 




46/38 



Table 5. Multiple Lorentzian fits to the optical PSDs from Night 1. See caption of Table O for details. Components marked with an asterisk are largely 
unconstrained. 
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